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CHAPTER  I 
INTRODUCTION 

A.  Background  Information 

In  the  design  of  a  marine  power  plant,  overall 
thermal  efficiency,  weight,  and  volume  are  the  driving 
concerns  (1] .  This  is  especially  true  for  the  steam 
plant,  which  is  the  predominant  form  of  propulsion  for 
large  naval  vessels.  While  a  large  marine  diesel  may 
provide  better  specific  fuel  consumption  than  a  fossil 
fired  steam  plant,  the  latter  retains  several  intrinsic 
characteristics  which  make  it  attractive  for  naval  pro¬ 
pulsion.  These  advantages  include: 

1.  Use  of  alternative  heat  sources. 

In  the  case  of  a  submarine,  aircraft  carrier, 
or  cruiser,  a  nuclear  reactor  provides  the 
thermal  energy.  In  smaller  ships,  a  boiler 
is  used  as  the  steam  generator. 

2.  Auxiliary  steam  generation. 

Auxiliary  steam  provides  water,  hotel  services 
and  drives  auxiliary  equipment  such  as  steam 
catapults  on  an  aircraft  carrier. 

3.  High  speed. 

The  steam  plant  has  a  much  wider  operating 
range  than  a  large  marine  diesel.  The  design 
of  a  naval  steam  plant  usually  orovides  for 
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redundancy  of  steam  generators  which  allows 
for  much  higher  speeds  than  a  comparable  sized 
diesel  plant. 

4 .  Endurance . 

The  steam  plant  which  has  a  nuclear  reactor 
as  its  heat  source  requires  neither  air  nor 
refueling  during  a  specific  mission. 

There  are  advantages  to  other  types  of  propulsion 
plants,  but  the  brief  list  above  points  out  that  for  cer¬ 
tain  applications  the  steam  plant  is  irreplaceable.  The 
high  speed,  large  shaft  horsepower,  and  endurance  required 
by  an  aircraft  carrier,  and  the  covert  endurance  required 
by  a  submarine  dictate  the  need  for  a  steam  propulsion 
system.  In  the  aircraft  carrier,  thermal  efficiency  may 
be  paramount  to  the  design  of  the  plant  due  to  large 
demands  by  auxiliary  systems,  but  in  the  submarine, 
weight  and  volume  are  critical. 

In  the  submarine  design,  25  per  cent  of  all  the  use¬ 
ful  weight  is  devoted  to  propulsion.  This  term,  for  the 
sum  of  all  useful  weight,  is  "normally  submerged  condi¬ 
tion,"  (NSC),  and  it  is  the  submerged  displacement  minus 
the  ballast  water  necessary  to  maintain  neutral  buoyancy 
[2] .  The  design  process  for  the  submarine  is  to  first 
lay  out  the  machinery  section  and  determine  its  critical 


length.  This  length  is  referred  to  as  the  "stack  length." 
Once  the  stack  length  has  been  established  along  with  its 
weight  and  volume,  the  remaining  required  weights  are 
added,  and  the  overall  buoyant  envelope  can  be  deter¬ 
mined  [3].  This  envelope  is  the  volume  necessary  for 
neutral  buoyancy.  The  stack  length,  therefore  seriously 
impacts  on  the  total  volume  of  the  submarine,  and  the 
weight  of  the  machinery  also  becomes  a  constraint.  With 
these  two  constraints  merely  removing  weight  from  within 
the  buoyant  envelope,  may  be  insufficient  to  allow  the 
reduction  of  the  envelope.  Weight  removal  must  often  be 
countered  with  addition  of  lead  ballast.  Conversely, 
reducing  the  volume  of  the  hull,  but  not  the  weight  causes 
negative  buoyancy.  Also,  the  decrease  in  volume  of 
specific  items  within  the  buoyant  envelope  allows  greater 
freedom  in  arrangement,  but  unless  accompanied  by  a  com¬ 
mensurate  weight  reduction,  a  decrease  in  NSC  will  not 
be  realized. 

Consider  the  fact  that  a  Los  Angeles  class  attack 
submarine  (SSN-688)  costs  460  million  FY-81  dollars  [4] . 
This  cost  includes  the  cost  of  construction  plus  all 
government  furnished  equipment.  The  dimensions  of  this 
submarine  are:  length,  362  ft;  diameter,  33  ft;  and 
submerged  displacement,  6900  tons  [  5] .  A  very  crude 


approximation  would  say  that  each  foot  of  length  costs 
one  million  dollars  and  displaces  19  tons.  So  a  moderate 
reduction  in  the  stack  length  which  is  accompanied  by  the 
appropriate  weight  reduction  should  be  worth  several  mil¬ 
lion  dollars  in  acquisition,  with  second  order  effects 
being  felt  in  horsepower  requirements. 


B.  Objectives  of  this  Work 

The  objective  of  this  thesis  is  to  develop  a  design 
methodology  for  an  alternative  surface  condenser  than 
those  presently  used  in  marine  steam  plants.  This  alter¬ 
native  condenser  is  expected  to  be  smaller  and  lighter 
than  existing  units,  and  it  will  have  a  higher  overall 
heat  transfer  coefficient. 

The  average  condensation  heat  transfer  coefficient 
on  a  smooth  surface  horizontal  tube  steam  condenser  is 
approximately  2000  BTU/hr  ft  °F.  Condensation  inundation 
of  tubes  which  are  lower  in  the  condensing  tube  bundle 
significantly  reduces  their  heat  transfer  effectiveness 
[6].  Vertical  orientation  of  the  smooth  tube  bundle 
would  not  solve  this  inundation  problem  because  a  developing 
boundary  layer  would  cause  similar  adverse  effects. 

Figure  1-1  shows  the  blocking  effect  of  the  fluid  for 
horizontal  and  vertical  tubes. 


In  1954  Gregorig  proposed  condensing  on  vertically 

oriented  tubes  with  fluted  axial  surfaces  [7] .  Surface 

tension  of  the  condensate  draws  the  liquid  into  the  valleys 

of  the  flute  where  it  then  drains  down,  leaving  a  thinner 

liquid  layer  on  the  peaks.  Figure  1-2  demonstrates  this 

concept.  Heat  transfer  coefficients  predicted  for  these 

condensing  surfaces  using  Gregorig' s  fluid  properties  and 

2 

temperature  drops  approach  8000  BTU/hr  ft  °F.  In  a  pro¬ 
cedure  developed  by  Zener  and  Lavi  [3]  and  refined  by 
Webb  [9]  for  the  design  of  an  optimized  Gregorig  conden¬ 
sing  surface,  heat  transfer  coefficients  as  high  as 
2 

36000  BTU/hr  ft  °F  are  reported  by  experiment  and  as  high 

“7 

as  55000  BTU/hr  ft  *F  are  predicted.  These  surfaces, 
however,  may  be  difficult  to  manufacture,  and  for  a 
generally  curved  repeating  surface  of  such  small  dimen¬ 
sions  as  are  required,  the  exact  geometry  may  not  be  a 
significant  influence  on  the  condensing  heat  transfer  co¬ 
efficient  [10] .  The  extrusion  capabilities  of  industry 
are  likely  to  set  the  limits  on  the  geometry  of  the  flute, 
and  for  these  reasons  a  quite  regular  geometry  such  as  a 
sine  flute  is  considered.  Longitudinal  machining  or  axial 
welding  to  achieve  a  specific  surface  geometry  is  not 
acceptable  for  condenser  tubes  which  are  used  in  naval 
vessels  [11] .  By  requiring  seamless  drawn  tubes  for 
naval  condensers,  residual  stresses  and  machining  flaws 


are  avoided  which  could  be  sources  of  unexpected  failure 
when  the  vessel  is  subjected  to  extremes  of  pressure, 
mechanical  shock,  or  thermal  shock. 

The  fluted  surface  having  a  sine  geometry  was 
analyzed  by  Yamamoto  and  Ishibachi  using  the  Gregorig 
model  [12] .  Their  analysis  enables  the  calculation  of  an 
average  condensate  film  thickness  on  the  surface  of  the 
flute.  The  problem  then  reduces  to  one  of  balancing  heat 
transfer  and  condensate  mass  flow.  The  phenomenon  of 
the  condensate  filling  the  valley  region  of  the  flute  is 
referred  to  as  "flooding,"  and  this  concern  for  flooding 
presents  the  same  problems  as  inundation  of  horizontal 
or  vertical  smooth  tubes.  Panchal  and  Bell  considered 
the  downstream  effects  of  condensation  as  the  valley  of 
the  flute  begins  to  fill  [13]. 

Proposing  an  alternative  condenser  design  which 
simply  re-orients  the  tubes  and  enhances  their  outer 
surfaces  should  accomplish  the  goals  of  being  smaller  and 
lighter.  Condenser  design  considering  alternative  mate¬ 
rials  must  be  treated  as  part  of  the  problem.  Materials 
which  warrant  consideration  are  Copper-Nickel  (Cu-Ni) 
90-10,  Cu-Ni  70-30,  and  Titanium,  commercially  pure, 
grade  2  (Ti,  CP,  Gr2) .  There  are  trade-offs  for  the  use 
of  these  materials,  the  most  important  of  which  are 
strength  and  weight  versus  thermal  conductivity. 


In  summary,  the  objective  of  this  thesis  is  to  pre¬ 
sent  a  methodology  for  a  marine  condenser  design  which  has 
vertical  fluted  tubes  and  is  possibly  fabricated  from  dif¬ 
ferent  materials.  This  condenser  must  be  lighter  and 
smaller  them  existing  condensers  to  warrant  its  manufacutr 
and  installation. 

C.  Design  Inputs 

Constraints  for  a  condenser  design  can  vary  widely. 
For  a  marine  condenser,  these  degrees  of  freedom  are 
somewhat  narrower,  and  the  design  process  is  governed 
by  the  concerns  for  weight,  volume,  and  thermal  efficiency 
For  naval  condensers,  reference  [11]  further  limits  the 
degrees  of  freedom  in  design.  These  limitations  arise 
because  of  the  requirements  for  the  naval  vessel  to  be 
able  to  operate  effectively  in  varying  environments. 

The  heat  transfer  requirements  are  for  the  conser¬ 
vation  of  energy  from  the  condensation  side  of  the  con¬ 
denser,  through  the  tube  wall,  and  into  the  cooling  water. 
Thus  the  heat  transfer  characteristics  can  be  analyzed 
separately  for  these  three  different  regions. 
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1.  Cooling  water. 

Various  correlations,  such  as  McAdams,  Colburn,  or 
Seider-Tate  can  be  used  to  determine  the  heat  transfer 
coefficient  for  turbulently  flowing  cooling  water  [6] . 
Fundamental  non-dimensional  groups  which  determine  the 
cooling  fluid  heat  transfer  coefficient  are  the  Reynolds 
number  and  the  Prandtl  number.  The  design  paprameters 
which  govern  these  numbers  are  temperature  dependent  fluid 
properties,  velocity,  mass  flow,  and  tube  hydraulic 
diameter. 

2.  Tube  Wall. 

Heat  transfer  resistance  through  the  tube  wall  is 
dependent  upon  tube  material,  wall  thickness,  and  a  scaling 
allowance. 

3.  Condensing  Fluid. 

The  heat  transfer  coefficient  on  the  condensing 
side  of  the  tube  will  be  shown  to  be  a  function  of  three 
non-dimensional  groups.  These  groups  are  specified  by 
surface  geometry,  condensing  pressure  (or  temperature) , 
and  condenstate  drainage  requirements  [7,  12,  13] . 

Structural  limitations  are  also  defined  in  reference 
[11] .  These  are  the  source  of  considerable  conflict  with 
the  heat  transfer  requirements.  The  desireable  high 
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strength  materials  generally  possess  lower  thermal  con¬ 
ductivities  and/or  inferior  metallurgical  performance  in 
salt  water  than  lower  strength  pure  alloys.  Table  1-1 
gives  a  brief  comparison  for  several  common  copper, 
nickel,  and  titanium  based  alloys.  Stiffners  are  also 
necessary  to  prevent  excessive  tube  vibrations.  These 
stiffners  may  be  used  to  provide  condensate  drainage 
strippers,  but  their  spacing  affects  the  condensation  heat 
transfer  coefficient  [13,  14] . 

Power  available  for  cooling  water  is  a  constraint 
which  also  affects  the  condenser  design.  Given  a  specific 
pump/motor  limitation,  by  either  size  or  power  available, 
the  velocity  and  mass  flow  of  coolant  are  constrained. 

Heat  losses  through  the  condenser  are  determined  by 
physical  shape  of  the  condenser,  i.e.,  number  of  passes, 
baffles,  etc.,  tube  length,  and  the  Reynolds  number.  As 
the  tube  length  of  the  condenser  decreases,  a  commensurate 
amount  of  pumping  power  becomes  available  for  either  in¬ 
creased  coolant  mass  flow,  or  increased  coolant  velocity 
[15]  . 

An  itemized  listing  of  design  parameters  for  a 
marine  condenser  can  be  specified  as  follows  [16] : 
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Table  1-1 

Thermal  Conductivity  and  Yield  Strength  of 
Several  Common  Metals  and  Select  Alloys 

[17,13] 


Metal/Alloy* 

Thermal 

Conductivity 

Yield 

Stress 

[BTU/hr  ft  ®F] 

[  ksi] 

Cu 

230 

*  4.8 

Cu-NI  90-10 

23 

16 

Cu-Ni  80-20 

21 

16 

Cu-Ni  70-30 

17 

25 

Ni  (98.9%  Pure) 

*** 

47 

16 

Monel  400  (66.5 

Ni-31.5 

Cu)  14 

45-25 

Monel  500  (66.5  Ni-29.5 
2. 7A1-0 . 6Ti) 

Cu-  10 

111 

Ti  (Cp-Gr  2) 

11.4-9.5 

40 

Ti-5Al-2.5Sn 

4. 5-4. 3 

120 

Ti-6A1-4V 

3.8 

120 

♦Properties  are  for  50°F<T<100°F,  annealed  condition. 

♦♦Yield  stress  is  defined  for  .2%  deformation. 

♦♦♦Nickel  and  Ni-Cu  alloys  are  unsuitable  for  marine  tubing 
because  they  are  subject  to  deep  pitting  in  seawater  at 
low  velocities  [23]  . 
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1.  Tube  material  and  material  properties. 

2.  Tube  geometry. 

a.  Inside  diameter,  outside  diameter,  and  wall 
thickness; 

b.  Flute  pitch  and  amplitude  and  nominal 
diameter . 

3.  Coolant  characteristics  . 

a.  Inlet  temperatures; 

b.  Mass  flow; 

c.  Velocity; 

d.  Physical  and  thermal  properties. 

4.  Condenser  operating  conditions. 

a.  Condenser  pressure; 

b.  Condensate  temperature  and  mass  flow; 

c.  Heat  transfer  requirements. 

5.  External  requirements. 

a.  Tube  length; 

b.  Tube  bundle  lay-out; 

c.  Number  of  passes; 

d.  Operating  depth  and  safety  factors; 

e.  Tube  sheet  and  stiffner  connections; 

f.  Air  ejector  locations; 

g.  Pumping  power  available  or  allowable  head 
loss. 


CHAPTER  XI 


TUBE  SELECTION 


A.  Materials 

The  materials  which  will  be  considered  for  marine 
condenser  tubing  are  the  copper-nickel  alloys,  Cu-Ni  90- 
10  and  Cu-Ni  70-30,  and  commercially  pure  titanium,  grade 
2  (Ti,  Cp) .  These  materials  are  of  interest  because  the 
U.S.  Navy  uses  both  Cu-Ni  alloys  [11] ,  and  Ti  is  a  possible 
alternative  material  which  will  provide  considerable 
weight  savings.  Table  II-l  lists  the  mechanical  and  ther¬ 
mal  properties  of  these  materials  along  with  their  compo¬ 
sitions  [17,18]. 

In  selecting  a  material  for  use  in  a  condenser, 
several  environmental  factors  must  be  considered.  These 
include  corrosion  characteristics  and  wear  properties  of 
the  material  in  a  salt  water  environment,  effects  of  bio¬ 
logical  fouling,  and  galvanic  coupling. 

1.  Cooling  water  velocity  effects. 

The  effects  of  sea  water  velocity  on  various  metals 
is  summarized  in  Figure  II-l. 

a.  Copper-Nickel  Alloys. 

Cooling  water  velocity  is  one  of  the  more  important 
factors  that  affect  Cu-Ni  condenser  tube  deterioration. 


Table  II- 


w  “  ' 

7T7« 

V'  l  - " 

•  «'■  ,r 

•’A  \ 

-  , 

'■ 

-23- 

A 

h 

iH  fH 

o 

m  oo 

<0  a 

1 

\ 

•  • 

■u  -u 

3 

VO 

o.  os  *r 

0  0 

1 

1 

-u  u 

o 

iH 

< 

to 

'—f 

o 

Ui 

X  X 

m 

aj 

<0  (0 

*» 

&4 

si 

s  e 

0 

• 

m  m  -h 

•u 

in  m 

X 

■u 

fH 

0  0  6* 

o 

•  • 

(0 

MM 

rH 

o  o 

o  o 

£ 

w 

m  (-•  i 

in  m  » 

ol 

m 

CN  fH  in 

H  H  (N 

CM 

\ 

• 

1  1 

• 

o 

OS 

C 

6*  6*  © 

X  X 

o 

6* 

0 

I  1  T3 

i0  It 

X 

-H 

*J  i4  i0 

Cl 

I  £  £ 

■u 

M  M  U 

z! 

1  o  o 

mm 

10 

Z  Z  C5 

•  • 

X 

h 

00 

c 

rH  fH 

Ifl 

0 

H 

O' 

^  ^  % 

£ 

£ 

•h 

so  vo  in 

©I 

o 

1 

X> 

in 

r- 

05 

O  >“t  VO 

hi 

co 

flj 

fH 

os  OS  CM 

H 

(U 

P*  CN 

X  X 

• 

Ul 

\ 

O  o  iH 

a 

i  i  l 

(0  10 

o 

o 

•  •  • 

a  o  oo 

Cl 

1  £  £ 

6* 

0) 

N 

1  o  o 

X 

fH 

z  z  z 

•  • 

<0 

e*  e*  6* 

•  fH  H 

X 

•h 

CO  C/3  C/3 

c 

(0 

CO 

u 

<  <  < 

0 

£ 

CO 

.  c 

cn  cn  m 

© 

-H 

m 

<0 

*H 

CN  CM  VO 

■U 

u 

C0  P» 

XI 

fH 

•H 

hi 

PS  PS  r* 

10 

•H 

•  • 

o 

u 

£ 

•  •  • 

z 

0) 

© 

1  ^  O 

• 

V* 

<•> 

0 

h< 

1  1  1 

o 

ai 

»H 

Jri 

Of 

.H  ** 

a. 

© 

s1 

• 

0 

a 

o 

o 

u 

c 

u 

X 

h 

© 

10 

•H 

■o 

£ 

cn 

CO 

vo  in  o 

e 

X  X 

Ul 

o 

5* 

X 

— i  cn 

0 

10  to 

•H 

u 

JQI 

1  £  £ 

• 

hi 

1  m  in 

o 

V 

o  o 

A 

a 

•  • 

*H 

•H 

o  o 

m 

M 

X 

1 

a 

in  ps  o 

(0 

<0 

u 

cn 

•  •  • 

z 

£ 

1 

+ 

P*  00  VO 

fH  PI 

CO 

o 

ps 

Zl 

o 

H 

•Hi 

1  1  1 

• 

Zl 

i  OS  OS 

o 

CM 

•H 

to 

i 

Stf 

a 

ps 

o  cn  in 

fH 

o  o 

£  £ 

•Hi 

£ 

1 

+ 

CM  CM  fH 

(0 

h  m 

3I 

©  (U 

6*1 

<u 

o 

4J 

i  i 

ul 

1  u  u 

1* 

fH 

© 

O  Q 

z 

os  p* 

\ 

> 

•h  -h  a* 

o  o 

o  o 

0 

z  z  u 

m  m 

H  cn 

•H 

1  1 

i  i 

i  i 

•H 

3  3  -H 

o  o 

o  os 

< 

U  U  6* 

OS  p* 

os  r- 

•H  «H 

a, 

•H  -H  h 

z  z 

u 

z  z  u 

3  3 

•H 

3  3  -H 

U  U 

6* 

U  U  6* 

atantt 


For  these  materials,  a  minimum  velocity  exists  below  which 
stagnation  occurs  and  crevice  corrosion  is  likely.  Above 
a  maximum  critical  velocity  pitting  and  erosion  occur. 

For  Cu-Ni  70-30  this  lower  velocity  limit  is  approximately 
2-3  ft/sec,  allowing  for  crevices  on  the  order  of  1.0  mil 
per  year  (mpy)  [19,20,21].  Crevice  corrosion  character¬ 
istics  for  Cu-Ni  90-10  are  similar  to  Cu-Ni  70-30  [22] . 
Various  values  are  reported  for  maximum  critical  veloci¬ 
ties,  above  which  impingement  attack  and  erosion  occur. 

The  maximum  sea  water  velocity  of  11  ft/sec,  specified 
by  reference  [11]  is  less  than  the  generally  accepted 
critical  velocity  for  the  Cu-Ni  alloys  which  is  12-15 
ft/ sec  [20,22]. 

b.  Titanium 

Ti  is  very  immune  to  the  effects  of  velocity.  This 
behavior  is  attributed  to  a  chemically  stable,  rapidly 
formed,  oxide  film.  Passivation  occurs  rapidly,  even 
in  the  presence  of  mechanical  damage  such  as  a  scratch. 
The  oxide  film  "repairs  itself"  virtually  instananeously 
[21,22,23].  In  sea  water  at  velocities  between  0-3  ft/ 
sec,  reported  general  corrosion  rates  are  nil  [21,23]. 
Crevice  corrosion  occurs  with  Ti  only  at  temperatures 
above  250°F  [23] ,  and  no  corrosion  occurs  at  sea  water 
velocities  up  to  25  ft/sec  and  possibly  greater  [21] . 
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2.  Biological  Fouling 

In  general,  the  more  corrosion  resistant  metals  and 
alloys  are  likely  to  experience  biological  fouling  under 
slow  moving  conditions.  Fouling  by  macro-organisms 
(mollusks,  barnacles,  etc.)  has  the  obvious  effect  of 
reducing  water  flow,  changing  velocity  profiles,  and  in¬ 
creasing  pumping  power  requirements.  Fouling  by  micro¬ 
organisms  (algae,  etc.)  has  the  effect  of  adding  another 
layer  of  thermal  resistance  to  the  inside  of  the  tube 
wall  [24] . 


a.  Copper-Nickel  Alloys 

Water  velocities  in  excess  of  about  5  ft/sec, 
temperatures  above  120 °F  and  below  50°F,  low  oxygen  con¬ 
tent  or  low  food  supply  (such  as  in  stagnant  or  closed 
systems)  in  the  water  will  effectively  prevent  fouling 
of  copper  alloys  by  macro-organisms  [20]  .  The  copper 
alloys  most  resistant  to  all  bio-foulding,  however,  are 
those  containing  more  than  85  percent  copper.  Thus,  even 
in  a  high  velocity  flow,  some  micro-biological  fouling 
(referred  to  as  "slime")  will  occur  on  Cu-Ni  70-30  sur¬ 
faces  due  to  bacterial  and  algae  growth  and  aaherance 

within  the  boundary  layer  [20,25,26,27].  The  irritant 

2+ 

or  toxic  effect  of  Cu  ions  is  believed  to  be  the  cause 
of  reduced  biological  fouling  on  Cu-Ni  alloy  tubes  [28] . 


b.  Titanium 


Since  Ti  is  such  an  immune  metal,  biological  fouling 
is  a  serious  problem.  The  threshold  velocity  to  prevent  . 
macro-biological  fouling  on  Ti  tubes  is  approximately  5 
ft/sec  [23,29].  Design  of  a  system  for  active  production 
of  toxic  agents  must  ensure  that  these  agents  are  not  en¬ 
vironmentally  hazardous.  A  chlorination  system  is  accept¬ 
able  to  the  environment,  but  the  production  of  chlorine 
ions  (Cl”)  is  accompanied  by  hydrogen  ions  (H+) .  Hydrogen 
ions  are  a  hazard  to  titanium  alloys,  and  except  In  the 
case  of  very  pure  metal  can  cause  stress  corrosion 
cracking  [30] . 

3 .  Galvanic  Coupling 

When  considering  condenser  materials,  in  the  presence 
of  a  strong  flowing  electrolyte,  namely  sea  water,  galvanic 
effects  are  significant.  The  material  itself,  with  its 
propensity  to  passivate,  and  the  velocity  of  the  sea  water 
dictate  the  kinetics  of  corrosion  due  to  galvanic  coupling. 
Titanium  is  a  much  more  active  metal  than  copper,  but  it 
is  because  of  this  activity  that  it  rapidly  forms  on  oxide 
and  becomes  passive.  Hence  it  is  more  noble  than  copper, 
and  its  general  corrosion  rates  are  much  less.  Table  II-2 
shows  the  relative  nobility  of  various  metals  in  flowing 
sea  water,  and  it  implies  that  Cu-Ni  piping  (used  for 
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Table  II- 2 

Galvanic  Series  of  Metals  in  Flowing  Sea  Water  [33] 


Anodic  or  Least  Noble 


Magnesium  and  magnesium  alloys 
Zinc 

Galvanized  steel 
Aluminum 
Cadmium 
Mild  steel 
Wrought  iron 
Cast  iron 

Stainess  steel  304  (active),  316  (active) 

Lead 

Tin 

Naval  brass  (60%  copper,  39%  zinc,  1%  tin) 
Copper 

Red  brass  (35%  copper,  15%  zinc) 

Copper-Nickel  90-10 
Copper-Nickel  70-30 
Nickel 

INCONEL  alloy  600  (78%  nickel,  13.5%  chromium, 

6%  iron) 

Nickel  aluminum  bronze 

Silver 

Titanium 

Stainless  steel  304  (passive) 

INCONEL  alloy  625 
HASTELLOY  alloy  C 
MONEL  alloy  400 
Stainless  steel  316  (passive) 

INCOLOY  alloy  825 

Graphite 

Platinum 

Cathodic  or  Most  Noble 


most  sea  water  piping)  would  deteriorate  in  the  proximity 
and  in  contact  with  Ti  [23,31].  Reference  [11]  requires 
zincs  as  cathodic  protection  in  Cu-Ni  tubed  condensers 
since  Cu  is  noble  relative  to  steel.  The  steel  is  noble 
relative  to  the  zincs ,  and  they  therefore  corrode  sacri- 
ficially.  Such  sacrificial  anodes  would  definitely  be 
required  in  a  Ti  tubed  condenser  [32].  In  fact,  the 
deterioration  rate  for  such  anodes  might  dictate  an 
entire  change  in  sea  water  piping  systems  to  a  more  noble 
metal  which  would  then  exacerbate  the  biological  fouling 
problem.  Figure  II-2  demonstrates  this  concern  for 
cathodic  protection  requirements  for  a  Ti  tube  condenser. 

The  significant  variables  for  material  use  in 
condenser  tubing,  apart  from  strength  requirements,  are 
sea  water  velocity,  galvanic  effects,  and  susceptibility 
to  biological  fouling.  The  latter  two  of  these  are  also 
functions  of  sea  water  velocity  as  shown  in  Figure  II-l. 

B.  Strength  Requirements 

Reference  [11]  specifies  for  surface  ship  conden¬ 
sers  a  minimum  tube  wall  thickness  of  .049  in.  with 
Cu-Ni  90-10  alloy.  The  tube  outside  diameter  is  also 
specified  as  5/8  (.625)  in.  Appendix  A  shows  the 


calculations  required  to  determine  condenser  tube  wall 
thickness.  These  results  are  presented  as  Table  II-3. 

The  minimum  wall  thickness  allowed  in  practice  for  Cu-Ni 
70-30  is  0.049  in.  and  for  Ti  is  0.035  in.  [34].  Thus,  from 
Table  II-3.  replacing  Cu-Ni  90-10  with  Cu-Ni  70-30  Ti  is 
most  practical  strictly  from  a  strength  consideration  for 
a  surface  ship.  Other  considerations,  such  as  increased 
erosion  protection  in  the  face  of  a  high  velocity  may 
force  the  change  in  materials  as  shown  in  Figure  ll-l. 

The  question  of  tube  wall  thickness  for  a  sub¬ 
marine  condenser  is  complicated  by  security  requirements. 
These  cam  be  treated  academically  by  the  following  con¬ 
servative  assumptions: 

1.  The  design  depth  of  a  submarine  is  not  the 
depth  at  which  it  operates.  A  safety  factor 
has  been  applied  which  may  be  either  additive 
or  multiplicative. 

2.  The  yield  criterion  for  tube  wall  material  also 
incorporates  a  factor  of  safety. 

Appendix  A  shows  the  calculations  required  to  pro¬ 
duce  Table  II-3.  Table  II-3  presents  various  condenser 
tubes  required  for  different  design  depths. 
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*  Table  n-3 

Tube  Specifications  for  Various  Design  Depths* 


Design  Depth  (ft) 

Material 

OD (in) 

t  (in) 

Gauge 

1050 

Cu-Ni  90-10 

.625 

.049 

18 

1050 

Cu-Ni  70-30 

.625 

.028 

22 

1050 

Ti  CP 

.625 

.019 

— 

1500 

Cu-Ni  90-10 

.625 

.083 

14 

1500 

Cu-Ni  70-30 

.625 

.049 

18 

1500 

Ti  CP 

.625 

.028 

22 

2000 

Cu-Ni  90-10 

.625 

.109 

12 

2000 

Cu-Ni  70-30 

.625 

.065 

16 

2000 

Ti  CP 

.625 

.035 

20 

2500 

Cu-Ni  70-30 

.625 

.083 

14 

2500 

Ti  CP 

.625 

.047 

18 

3000 

Cu-Ni  70-30 

.625 

.095 

13+ 

3000 

Ti  CP 

.625 

.058 

17+ 

1050 

Cu-Ni  90-10 

.500 

.035 

20 

1050 

Cu-Ni  70-30 

.500 

.028 

22 

1050 

Ti  CP 

.500 

— 

— 

1500 

Cu-Ni  90-10 

.500 

.065 

16 

1500 

Cu-Ni  70-30 

.500 

.035 

20 

1500 

Ti  CP 

.500 

.022 

— 

2000 

Cu-Ni  90-10 

.500 

.083 

14 

2000 

Cu-Ni  70-30 

.500 

.049 

18 

2000 

Ti-  CP 

.500 

.028 

22 

2500 

Cu-Ni  70-30 

.500 

.065 

16 

2500 

Ti  CP 

.500 

.035 

20 

3000 

Cu-Ni  70-30 

.500 

.083 

14 

3000 

Ti  CP 

.500 

.049 

18 

Thickness  for  tubes  has  been  rounded  to  the  nearest 
even  gauge  except  as  noted  by  +  . 


+  + 
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CHAPTER  III 

* 

ELEMENTS  OF  HEAT  TRANSFER 


Nomenclature 

a  amplitude  of  the  flute  (ft) 

A  area  (ft2) 

B  dimensional  group  defined  in  eq.  (III-7) 

(f t/°F) 

C,  Cp  specific  heat  (BTU/lbm°F) 

D  diameter  (ft) 

F  number  of  flutes  on  the  tube  surface 

2 

g  gravitational  constant  (lbm  ft/lbf  hr  ) 

c 

G  mass  flux  (lbm/ft2  hr) 

h  heat  transfer  coefficient  (BTU/hr  ft2°F) 

hjg  latent  heat  (BTU/lbra) 

K  thermal  conductivity  (BTU/hr  ft°F) 

L  length  (ft) 

M  non-dimensional  group  defined  in  eq.  (III-16) 

Nu  Nusselt  number 

p  pitch  of  the  flute  (ft) 

2 

P  pressure  (lfb/ft  )  or  (in-hg-abs) 

Pr  Prandtl  number 

q  heat  flow  (BTU/hr) 

r  radius  of  curvature  of  the  fluted  surface 


(ft) 


-34- 


R  thermal  resistance  (hr  ft2°F/BTU) 

Re  Reynolds  number 

Sv  flute  half  perimeter  (ft) 

Sc  flute  arc  length  on  which  condensing  occurs 

(ft) 

t  tube  wall  thickness  (ft) 

T  temperature  (°F) 

AT  Ts  -  Tw  (°F) 

0  overall  heat  transfer  coefficient  (BTU/ 

hr  ft2°F) 

V  velocity  (ft/hr) 

W  axial  mass  flow  of  liquid  (lbm/hr) 

s 

y  ■  tube  surface  coordinates 

z 

4 

Greek  letters 

a  height  of  the  condensate  in  the  center  of 

the  flute  (ft) 

r  mass  flow  rate  in  S  direction  per  unit 

length  (lbm/hr  ft) 

<5  condensate  film  thickness  (ft) 

angular  representation  of  a  point  along 
the  flute  surface 


8 
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y  dynamic  viscosity  (lbm/hr  ft) 

p  density  (lbm/ft3) 

a  surface  tension  (Ibf/ft) 

4  non-dimensional  group  defined  in  eq.  (B-12) 

ft  non-dimensional  group  defined  in  eq.  (III-7) 

Subscripts 

b  at  bulk  temperature 

c  coolant 

D  referenced  to  hydraulic  diameter 
f  at  flooding  conditions 

i  inside 

n  nominal 

o  outside 

s  saturation 

w  wall 

sc  scale 

Symbols 


average 


A.  Heat  Transfer  Coefficient  of  Cooling  Water 

No  simple  analytical  solution  exists  for  heat  trans¬ 
fer  in  turbulent  pipe  flow.  The  McAdams  correlation  (or 
Dittus-Boelter  equation)  is  widely  accepted  for  the 
determination  of  the  heat  transfer  coefficient  for  the 
cooling  water  [6,39]. 


»  0.023  pleDj0*8(Pr)£  (Ill-la) 

where 

1)  n  *  0.3  if  the  fluid  is  being  cooled 

*  0.4  if  the  fluid  is  being  heated 

2)  all  fluid  properties  are  evaluated  at  the 
average  (or  bulk)  fluid  temperature 

3)  2300  <  RBq  <  10^  where 

ReD  =  Reynolds  number  based  upon  the 
hydraulic  diameter 

4)  0.5  <  Pr^  <  120  where 

Pr^  =  Prandtl  number  based  on  bulk  temperature 

5)  L/D  >  50 
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For  flow  through  a  condenser  tube,  these  conditions 
are  usually  met  and  equation  (III-l)  is  used  as  a  basic 

i 

relationship. 


B.  Heat  Transfer  Coefficient  for  a  Tube  Wall 

For  materials  such  as  pure  copper  which  have 
extremely  high  values  for  thermal  conductivity,  the 
thermal  resistance  of  the  condenser  tube  wall  is  almost 
negligible.  From  Table  1-1,  it  is  apparent  that  the 
strength  requirements  for  marine  condenser  tubes  mandate 
the  use  of  materials  with  relatively  low  thermal  con¬ 
ductivities  . 

The  thermal  resistance  of  the  tube  wall  is  the 
reciprocal  of  the  heat  transfer  coefficient  of  the  wall. 
For  a  smooth  (non-enhanced)  tube,  this  may  be  written  as: 


x  zn(Do/Di)D 

K  *  1 


(IXI-2) 


For  a  fluted  tube,  however,  the  nominal  diameter  is  used 
in  place  of  the  outside  diameter  [40] .  The  nominal 
diameter  is  defined  as: 


D_  a  D  +  2a 
n  o 


( III-3 ) 


This  is  shown  in  Figure  Ill-l.  Use  of  the  nominal  dia¬ 
meter  will  be  shown  to  be  reasonable  for  two  reasons: 

1.  The  size  of  the  flute  is  small  with  respect 
to  wall  thickness .  (This  is  exaggerated  in 
Figure  III-l) . 

2.  Heat  transfer  is  blocked  by  condensate  at  the 
bottom  of  the  flute  valley  as  shown  in  Figure 
III-2.  [91. 

Furthermore,  it  is  convenient  to  work  with  the  nominal 
diameter  vice  the  outside  perimeter  when  dealing  with 
the  heat  flux  and  condensing  heat  transfer  coefficient. 
Equation  (III-2)  can  be  written  as: 

-  \  - -  (IH-2a) 

*n(VDi>  Dn 

K  2 


C.  Heat  Transfer  Coefficient  for  Condensation 

The  coordinate  system  for  this  analysis  is  schemati¬ 
cally  depicted  in  Figure  III-2.  The  arc  length,  S,  is 
calculated  by 


f  1  (1+  fill] 

)  2ir  r  [  p  j 


.  2  ) 
s  m  s  J 


1/2 


dd> 


(III-3) 


S  - 


[• 


=  u  s  r  r  &e  -  6  si  a  e/2  x  2 

3  0  =  6  co 5 e/i  “  6  cos  e/a 


,  m-3-  S c. Vj  tm  3tvC  of  fluid  Sur  face 
j  force,  balance  between  pressure 


w  r  fac€.  ten  s  i  c  n  . 
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Table  III-l 

Sample  Calculations  for  Condensate  Film  Thickness 


a/p 

.15 

.35 

.45 

in 

in 

• 

6_X105 

n 

ft 

0 

1.5811 

.  4878 

.3262 

.2422 

5 

1.7219 

.6673 

.5436 

.4802 

10 

2.0708 

1.1451 

1.0717 

1.0522 

15 

2.5861 

1.8466 

1.8404 

1.8785 

20 

3.1842 

2.6712 

2.7420 

2.8451 

25 

3.7202 

3.4183 

3.5596 

3.7215 

30 

4.0099 

3.7606 

3.9258 

4.1083 

35 

3.9580 

3.5078 

3.6242 

3.7693 

40 

3.6563 

2.8574 

2.8858 

2.9616 

45 

3.2976 

2.1107 

2.0450 

2.0468 

50 

3.0909 

1.4664 

1.3067 

1.2374 

55 

3.4029 

1.1087 

.8361 

.6900 

60 

6.7751 

1.8444 

1.2311 

.8894 

a  -  2.00  x  10‘3  ft 

0  *  3.7649  x  10"9 

49  »  :r/60 

5  *  f(nA9) 

n 


From  the  results  of  equation  (III-6) ,  solved  for 


the  arc  length  Sc  ■  Sv,  NuQ  can  be  calculated  as  shown 

1/4 

in  Appendix  B.  A  non-dimensional  number,  NuqQ  ,  can 
be  calculated,  and  from  equation  (III-6)  and  the  defini¬ 
tion  of  Nu  , 
o 


NuQa1/4  *  f(a/p) 


{ III-8 ) 


This  result  is  presented  in  Table  III-2  and  plotted  in 
Figure  I1I-4. 

At  a  distance  Z  down  the  tube,  condensate  W  runs 

down  the  flute,  and  the  thickness  of  this  condensate,  a 

(Figure  III-2) ,  increases  with  Z.  Here  the  heat  transfer 

is  neglected  in  the  valley  of  liquid  thickness  a,  and  the 

integration  of  equation  (III-6)  stops  at  S  *  S  .  S  is 

c  c 

determined  by  equation  (III- 3)  and  9  . 

.  v» 


e 


c 


(III-9) 


When  a  »  2a  the  valley  is  flooded  [13] ,  and  w  *  W  _ 

where 


w 


f 


3.33 


(III-10) 
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Table  III-2 

1/4 

Sample  Calculations  for  Nu  a  =  f(a/o) 

o 


axlO3 

(ft) 

a/p 

Nu 

o 

axio9 

Nu  21//4 
O 

10.0 

.1 

532.31 

3.7649 

4.1697 

5.0 

.1 

532.31 

3.7649 

4.1697 

1.0 

.1 

532.31 

3.7649 

4.1697 

2.0 

.15 

553.65 

3.7649 

4.3363 

10.0 

.25 

326.54 

37.649 

4.5486 

10.0 

.25 

580.68 

3.7649 

4.5436 

10.0 

.25 

1032.62 

.37649 

4.5486 

2.0 

.25 

580.69 

3.7649 

4.5486 

2.0 

.40 

595.41 

3.7649 

4.6639 

2.0 

.45 

596.89 

3.7649 

4.6755 

2.0 

.50 

597.45 

3.7649 

4.6799 

2.0 

.55 

597.36 

3.7649 

4.6792 

2.0 

.60 

596.80 

3.7649 

4.6749 

2.0 

.65 

595.92 

3.7649 

4.6680 

2.0 

.70 

594.81 

3.7649 

4.6593 

2.0 

.75 

593.55 

3.7649 

4.6494 

or 


wf  *  f  9  V  (a)3  p  exp 


-3.33 


r  > 

a 

Ip. 


(Ill.lOa) 


A  Reynolds  number  is  defined  [13]  as 


Re 


4  W 

5  u 
v 


(III. 11) 


Calculations  by  Panchal  and  Bell  [13]  show 


Wj  Re^ 


<3 

2a 


3.252 


(III. 12) 


From  Appendix  B  and  equations  (III-9)  and  (III-12) 


h*  >  Nuz  .  , 

a 

Wz 

a 

ho  Nuo 

Re- 

'  p 

j 

» f 

Wf  ' 

p 

(III. 13) 


This  is  shown  by  example  in  Figure  III-5. 

An  energy  balance  on  the  tube  surface  requires 


tD_ 


h_dZ  it  D.  AT  ■  h-  dW 
z  n  fg 


n 


(III. 14) 


This  equation  is  re-arranged  and  integrated  between  Z  *  0 
and  Z  »  L. 


L  p  AT 


V 

K 


(III. 15) 


K  Nu  a1/4 
-E. _ o _ 


hfg  Wf  NUZ 


BAT 

a 


1/4 


Nu 

o 


dW 

Nuz 


By  grouping  the  constants  in  the  right-hand  expression, 
M  can  be  defined  as 


K(a) 


1/4 


M  5 


Nu  n 
o 


1/4 


hfg  (B) 


ITT 


LAT 


3/4 


W, 


Figure  III-5  is  replotted  as 


(W/Wf)L 

0 


(III. 16) 


Nu 

O 

Nuz 


(III. 17) 


and  presented  as  Figure  III. 6.  The  integral  in 
equation  (III. 16)  is  evaluated  from  Figure  III-6,  and 
two  of  the  cases  are  presented  in  Figures  III-7a  and 
III-7b.  These  results  are  also  tabulated  and  presented 
as  Table  III-3.  Note  that  the  integration  smooths  out 
most  of  the  dependency  on  (a/p)  so  that 


0  10 


20  30  40  50  60  70  80  90  100 

W/Wf  X  10 


Figure  III-7b:  M  vs.  W/Wf 
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Table  III-3 


w 

Calculated  Values  for  M  and 

f 

for  Several  Different  Values  of  a/p 


M  x  103 


x  IQ' 


a/p  ■  .15 


a/p  =  .55 


— l 


0 

0 

1 

1.06 

2 

2.18 

3 

3.35 

4 

4.44 

5 

5.61 

6 

6.79 

7 

7.98 

8 

9.19 

9 

10.40 

10 

11.63 

15 

17.87 

20 

24.36 

30 

37.61 

40 

51.48 

50 

65.97 

60 

80.75 

70 

96.07 

30 

111.63 

90 

127.67 

100 

143.96 

110 

160.70 

0 

1.07 

2.21 

3.38 

4.56 

5.77 

6.97 

8.22 

9.46 

10.72 

11.98 

18.43 

25.04 

37.75 

52.56 

67.14 
81.95 

97.15 
112.55 
128.30 
144.21 
160.46 


0 

1.07 

2.24 

3.42 

4.62 

5.85 

7.08 

8.32 

9.58 

10.85 

12.12 

18.62 

25.26 

38.97 

52.76 

67.20 

81.88 

96.90 

112.08 

127.55 

143.16 

159.05 


M  »  f  (W/Wf) 


(III. 18) 


M  ji  f  (W/Wf/  a/p) 


(111.13a) 


The  data  in  Table  III-3  suggests  a  relationship  for  M 
such  that 


M  =  1.7312  (W/Wf) 


1.0839 


(III. 19) 


for  W/Wf  >_  0.003 


The  average  heat  transfer  coefficient,  K,  is 
defined  as 


U  .  = 


(q/an> 


AT 


(III. 20) 


which  may  be  written  as 


_  •>£„  »£  ("J 

h  *  pTlT 


h-  W.  .  6027  (M) 
_ £ 


.9226 


p  L  AT 


(III. 21) 


Substituting  for  M  from  equation  (III. 16)  gives 


where  f(a/p)  *  Nu  from  Figure  III-4.  At  L  =  0 

this  should  show  that  h  *  K  .  The  relationship  for  w /w^ 

o  - 

amd  M  fails  for  w/w^  <  0.003,  and  this  equation  for  h 
goes  to  «. 


D.  Application  to  Condenser  Design 

In  the  design  of  a  fluted  condenser,  it  may  be 
desireable  to  place  stripper  plates  along  the  length  to 
remove  the  condensate  well  before  it  reaches  the  flooding 
flow  rate  -  perhaps  keeping  w/wg  <_  0.1.  Figure  III-8 
shows  such  stripper  plates. 

To  determine  L  amd  w/Wg,  a  system  of  4  equations 
must  be  solved.  Other  unknowns  to  this  system  of  equa¬ 
tions  are:  (q/An) •  Tc2  and  Therefore,  one  of  these 

five  variables  must  be  specified.  In  practice,  either 
L  or  w/Wg  would  generally  be  chosen.  The  system  of 
equations  is: 


Figure  U[-S  -  Schematic  diagram  showm 
Stripper  plate  spacing  and  axial 
temperature  distri  bu  t  ion. 
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(q/v 


hf(?  (w/w£)  (AT) 
P [L  (AT)3/4] 


ran  3/4 


(III. 23) 


which  comes  from  the  conservation  of  energy  on  one  flute 
of  the  tube  surface. 


(q/An) 


Wc  Cc  (Tc2  "  Tcl> 
*  Dn  L 


(III. 24) 


which  comes  from  the  conservation  of  “nergy  in  the  cooling 
water . 


(q/An) 


(Tc2-Tc1) 
It  -at-tc,  f 

3  i 

'  (T  -AT-Tc, ) 

_  S  i  _ 


(III. 25) 


which  is  the  conservation  of  energy  through  the  tube  wall 
using  an  overall  heat  transfer  coefficient  and  log  mean 
temperature  difference  for  the  coolant. 


1  ,  n  1  ,  1 

Un,w  Di  hc  hsc 


in  (D/D. )  D 
n  i  n 


(III. 26) 
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For  short  lengths  or  small  temperature  rises  in  the  cool 
ant,  the  log  mean  temperature  difference  may  be  replaced 
by 


in 


(Tc2  -  Tc1) 

"(T^-AT-Tc^)- 

(T  -AT-Tc~) 

.  s  2  _i 


(III. 27) 


where  Tc  =  j  (Tc1  +  Tc2)  (III. 28) 

The  fourth  equation  of  the  system  comes  from  the  defini¬ 
tion  of  M  in  equation  (III. 16) 


4) 


M  wf (B)1/4  hf 
(NuQn1/4)  (a) 


K 


(III. 29) 


where  M  is  approximated  by  equation  (III. 19). 

With  L  specified  and  w/w^  allowed  to  float,  the 
solution  to  these  equations  is  as  follows: 

1)  Assume  Tc2 

2)  Assume  AT 

3)  Solve  for  (q/A^)  in  equation  (III. 25) 

4)  Solve  for  (w/wf)  in  equation  (III. 23) 

5)  Solve  for  M  in  equation  (III. 19) 


6) 


Solve  for  IT  in  equation  (III. 29) 
Repeat  steps  3-6  until 


7) 


ATi  ‘  ATi+l'  ± 


’AT 


8) 

9) 


Solve  for  TC2  in  equation  (III. 24) 
Repeat  steps  2-8  until 

lTc2j  -  Tc2j+1|  i  sTc  • 


With  w/w^  specified  and  L  allowed  to  float,  the 
solution  to  these  equations  is  as  follows: 

1)  Assume  Tc2 

2)  Assume  L 

3)  Solve  for  AT  in  equation  (III. 29) 

4)  Solve  for  (q/An)  in  equation  (III. 25) 

5)  Solve  for  L  in  equation  (III. 23) 

6)  Repeat  steps  3-5  until 

!Li  "  Li+l'  -  eL 


7)  Solve  for  Tc2  in  equation  (III. 24) 

8)  Repeat  steps  2-7  until 

I Tc2 j  "  Tc2j+ll  -  £Tc2’ 

Table  III-4  shows  the  results  of  these  calculations 
for  a  typical  condenser,  neglecting  thermal  resistance 
in  the  wall.  Appendix  C  shows  the  calculations  for  the 
solution  to  this  system  of  equations  as  applied  to  a 
condenser  design. 


Table  III-4 


Sample  Results  for  a  Fluted  Condenser* 

Given:  DQ  *  .625  in 

t  »  .049  in 

Tcx  -  75 . 0°F 

v  *  7.5  ft/sec 
c 

w  *  12000  GPM  (6.1513  x  106  Ibm/hr) 
c 

a  *  1.50  x  10~3  ft 

a/p  *  .31205 

P  *  2.5  in-Hg-abs 

G 

Nu  Q1/4  *  4.612 

o 

B  »  4.3777  x  10'12ft/°F 

D  *  5.5083  x  10“2  ft 

n 

P  «  4.8069  x  10“3  ft 

F  -  36 

No.  of  tubes  *  2354 

*  5.4020  lbm/hr 


Section 

1 

2 

3 

4 

AT(#F) 

1.084 

0.764 

0.523 

0.350 

w/wf 

0.10 

0.079 

0.060 

0.046 

L  (ft) 

3.252 

3.252 

3.252 

3.252 

(q/An) (BTU/hr f t2) 

35  650 

23000 

21550 

16300 

TCj 

83.17 

89.59 

94.53 

98.26 

«total  *  l*34s  x  108  BTU/hr 

♦wall  resistance  neglected;  cooling  fluid  is  sea  water. 
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E.  Heat  Transfer  Resistance  Due  to  Scaling 

Allowance  for  dirt  or  scale  can  be  made  by  either 
including  separate  thermal  resistances  for  the  scale  or 
fouling  associated  with  the  tube  material  and  fluid  or  bv 
allowing  a  percentage  reduction  of  the  overall  heat  trans¬ 
fer  coefficient,  U.  Reference  [11]  requires  the  latter 
method  using  a  15  percent  reduction.  Reference  [6]  recom¬ 
mends  a  straight  allowance,  regardless  of  material  but 
dependent  upon  cooling  medium.  The  value  recommended  for 
sea  water  is: 


.0005  (hr  ft2°F/BTU) 


Reference  [41]  recommends  including  a  separate  resistance 
for  metal  oxide  on  each  side  of  the  tube  and  a  separate 
fouling  resistance.  This  gives  a  value  for  U  of  the  form. 


_1 

h 


condensing 


D^’fT  +RW  + 

i  sea  water 


ox  +  R* 

(III-31) 


Values  for  Roxj_^e  are  listed  in  Table  III-5.  Since 
Rfouiing  i®  a  virtual  unknown  and  the  quantity  2Rqx  +  R* 
can  be  approximated  for  all  cases  in  sea  water  as 
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2R  +  R*  -  Rc_  "-.0005  (hr  ft  °F ) /3TU) 

OX  sc 


This  gives  a  value  for  Un  of  the  form 


_1_  ,  1  +  __n 

Un  ^condensing  Di 


in 


V 


i  + 1 

sea  water  scale 


n 


(III. 32} 


tube 


where  h  -  2000  BTU/hr  ft‘°F. 

SC 


TABLE  III-5 


Values  of  Rqx,  the  Resistance  of  the 
Oxide  Film  on  the  Inside  and  Outside 
Surfaces  of  Clean  Tubes  [41] 


Tube  Material 


oxide 


hrft  °F 


BTU 


i 


Admiralty  Metal 
Aluminum  Bronze 
Aluminum  Brass 
Cu-Ni  90-10 
Cu-Ni  80-20 
Cu-Ni  70-30 


0.000136 

0.000153 

0.000167 

0.000178 

0.000193 

0.000243 


Titanium 


0.000195 


CHAPTER  IV 


DESIGN  PROPOSALS 


A.  Condenser  Geometry 

A  primary  objective  of  this  work  is  to  provide  a 
methodology  for  a  condenser  design  which  allows  the  naval 
architect  to  specify  condenser  length.  The  naval  archi¬ 
tect  is  constrained  by  weight  and  volume,  but  the  dimen¬ 
sion  which  is  most  critical  in  marine  plant  layout  is  the 
tube  length.  The  transverse  area  of  the  tube  bundle  can 
grow  with  significantly  less  impact  on  the  other  en¬ 
gineering  systems  than  axial  growth.  By  specifying  the 
tube  length,  the  maximum  tube  deflections  and  vibration 
characteristics  can  be  predetermined.  This  allows  the 
placement  of  tube  and  shell  stiffners  which  may  also  act 
as  condensate  stripper  plates  as  suggested  in  Figure 
III-8.  For  typical  spacing  requirements  of  1-8  feet, 
flutes  with  amplitude  of  12-24  mils  (1.0  -  2.0  x  10“3  ft) 
are  predicted  by  this  work.  Appendix  C  demonstrates  how 
a  length  may  be  specified  for  a  condenser  design.  A 
typical  double-pass  fluted  condenser  drawing  is  shown 
in  Figure  IV-1. 


Figure -IV-]*:  Poaaible  configuration  for 
a  2  Pats  Vertical  Condenser 


B.  Tube  Attachment 

A  serious  consideration  for  a  naval  condenser  is 
the  ability  to  pull  and  replace  or  plug  tubes  without 
having  to  remove  the  entire  condenser  or  tube  sheet. 
Tubes  in  naval  condensers  are  usually  attached  to  the 
tube  sheets  by  rolling  their  ends  and  actually  pressing 
them  into  the  tube  sheet  [41,42].  This  procedure  is 
well  established  for  Cu-Ni  alloys,  and  rolled  Ti  tube- 
to-tube  sheet  joining  procedures  are  being  established 
to  meet  naval  heat  exchanger  requirements  [43] . 

Figure  IV-2  is  a  drawing  of  a  proposal  to  facili¬ 
tate  fluted  tube  attachment  to  the  tube  sheet  while 
maintaining  the  ability  to  pull  a  tube  if  required. 

As  shown  in  Figure  IV-2,  a  forged  collar  can  be  placed 
around  the  fluted  surface  at  points  where  attachment  is 
desired.  This  collar  could  be  snugged  into  place  by 
heating  it,  prior  to  attachment,  and  allowing  it  to  cool 
around  the  tube.  Reference  [41]  shows  a  similar  attach¬ 
ment  for  ferrules  to  allow  longitudinal  expansion. 
Reference  [42]  reports  successful  rolling  of  5/8  inch 
diameter  Ti  tubes  with  0.049  inch  wall  thickness.  The 
composite  wall  thickness  for  a  Ti  tube  with  flute  ampli¬ 
tude  of  18  mils  and  tube  wall  thickness  of  0.035  inches 
is  0.710  inches.  This  may  pose  a  problem,  and  further 
research  is  warranted. 


CHAPTER  V 


RESULTS  AND  CONCLUSIONS 


A.  Comparison  with  Nusselt  analysis 

The  average  heat  transfer  coefficient  for  conden¬ 
sation  on  a  fluted  surface  is  given  by  equation  (III. 22) 


"  3  .  “1.2307  .  .  .2307 

K  phf  (ogc) 


.6027 


.2037  ...  .  n.9226 

x  a  [f  (a/p)J 


hfqWf~|  • 
_  LAT  _ 


0774 


(III. 22) 


Prom  the  Nusselt  analysis  for  condensation  on  a 
smooth  vertical  surface. 


.943 


K3ph, 


(p-p  )g 


(V.l) 


These  two  equations  show  a  much  different  dependence 
on  the  axial  length  L,  but  h^  is  strongly  influenced  by 
pitch.  AT  also  affects  h^  more  than  hNu.  All  other 
things  being  of  a  similar  magnitude  or  extremely  weak 
functions,  is  dominated  by  p,  AT,  and  a.  The  Nusselt 

analysis  shows  that  h^  is  strongly  affected  by  L  and  AT. 


V  %*  V  ”  -  •’ 
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Figure  V-l  is  a  graph  relating  the  average  heat 

transfer  coefficient  with  length.  Since  M  as  defined  in 

equation  (III. 16)  is  linear  with  L,  the  relationship 

between  (w/w^)  ,  M,  and  L  for  a  given  h^  (or  hL/hQ  = 

Nu_/Nu  )  can  be  seen  to  be  uniaue. 

L  o 


B.  Conden'  Comparison  Using  Titanium 

To  cc  ..^are  a  Ti  condenser  with  the  Cu-Ni  70-30 
condenser  cited  in  Appendix  C,  another  constraint  must 
be  included.  Since  the  Ti  tubes  have  thinner  walls, 
less  of  them  will  be  required  for  the  sea  water  mass 
flow.  This  implies  that  the  Ti  condenser  will  either 
be  longer  or  transfer  less  heat.  An  alternative  is  to 
allow  the  mass  flow  of  the  coolant  to  increase,  i.e., 
increase  the  number  of  tubes.  This  will  be  the  approach 
taken  here  so  that  the  condenser  length  and  overall  heat 
transfer  remain  the  same.  Table  V-l  shows  the  differences 
between  these  two  condensers. 

No  attempt  has  been  made  in  Table  v-l  to  assess 
the  overall  weight  change  between  Cu-Ni  70-30  and  Ti  as 
condenser  materials  for  the  example  condenser.  This 
overall  weight  change  must  include  tube  sheets,  structural 
elements,  waterbox  heads,  affiliated  piping,  foundations, 
etc.  Assessment  of  these  weights  is  beyond  the  scope  of 
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Table  V-l 

Comparison  Between  Condensers  Using  Cu-Ni  and  Ti 


Cu-Ni  70-30 

Ti 

0o  (in) 

.625 

.625 

t  (in) 

.065 

.035 

a  (mils) 

18 

18 

a/p 

.13869 

.13869 

Tc±  (°F) 

65.0 

65.0 

V  (ft/sec) 

C 

8.0 

8.0 

Wc  (GPM) 

7900 

8600 

nt 

3294 

2852 

P_  (in-Hg-abs) 
c 

5.0 

5.0 

h  (BTU/hr  ft2°F) 

SC 

2000 

2000 

K  (BTU/hr  ft°F) 

17.0 

11.4 

F 

16 

16 

No.  of  passes 

2 

2 

tIV1/4 

4.305 

4.305 

Section  1 

Econ<J  (BTD/hr  ft2°F) 

12312 

11855 

Un  (BTU/hr  ft2°F) 

487.7 

535.7 

4T(°F) 

2.489 

2.895 

w/wf 

0.060 

0.060 

L  (ft) 

4.957 

4.426 
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Table  V-l  (Continued) 


Cu-Ni  70-30 


q/An  (BTU/hrft*) 


30649 

76.36 


34325 


74.03 


Section  2 


5cond  <BTU/hrft^F) 
On  (BTU/hrft2(*F) 

AT  ( °F) 


13197 

501.3 


1.988 


12506 

547.9 


2.435 


0.051 


0.053 


L  (ft) 


4.957 


4.426 


q/An  (BTU/hrft) 


Tc2  (°F) 


Section  3 


^cond  (BTU/hrft2 °F) 
Un  (BTU/hrft2°F) 

AT  (°F) 


26230 

86.08 


15302 

510.8 


1.528 


30451 

82.05 


13628 


557.6 


2.007 


L  (ft) 

q/An  (BTU/hrft2) 
Tc2  CF) 


Section  4 


Soond  <BTU/hrft2"F) 

U  (BTU/hrft2 ®F) 
n 

AT  (°F) 


0.019 

2.086 

23383 

89.72 


15800 

514.9 


1.377 


0.034 

3.148 

27351 

87.16 


14205 

564.2 


1.754 


0.018 


0.031 


-73 


Table  V-I  (Continued) 


Cu-Ni  70-30 

Ti 

L  (ft) 

2.086 

3.148 

q/An  (BTU/hrf t2) 

21759 

24919 

Tc2  (°F) 

93.12 

91.83 

Section  5 

Econd  <BTO/nuft2«F) 

15178 

14405 

Un  (BTU/hr  ft2 °F) 

519.0 

569.9 

AT  (°F) 

1.262 

1.539 

w/wf 

0.038 

0.039 

L  (ft) 

4.957 

4.426 

q/An  (BTU/hrf t2) 

19160 

22166 

Tc2  (°F) 

100.22 

97.66 

Section  6 

Eoond  (BTU/hrft2°P) 

16441 

15326 

Un  (BTU/hrft2°F) 

525.9 

576.7 

AT  (#F) 

0.974 

1.258 

w/wf 

0.031 

0.034 

L  (ft) 

4.957 

4.426 

q/An  (BTU/ft2) 

16011 

19287 

Tc2  (°F) 

106.15 

102.7 

«total  (BTD/hr) 

1.57  x  108 

1.57  x  10 

Bundle  Diameter  (in) 

61.8 

57.0 

Tube  Bundle  Weight  (tons) 

10.31 

3.78 

Weight  of  Added  Water  (tons) 

— 

2.67 
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this  work.  Table  V-I  shows  that  a  smaller  lighter  con¬ 
denser  can  be  designed  using  Ti ,  but  the  penalty  is  a 
requirement  for  increased  coolant  flow.  Table  v-2  is 
an  estimate  of  a  comparable  horizontal  condenser  using 
Cu-Ni  70-30  and  designed  to  meet  the  same  requirements 
as  the  condensers  presented  in  Table  V-l. 

C.  Conclusions 

This  thesis  has  presented  a  design  methodology  for 
a  fluted  condenser.  It  has  shown  by  example  in  Table 
V-2  that  a  fluted  condenser  can  be  designed  which  is 
more  compact  than  a  horizontal  condenser.  Overall  den¬ 
sity  of  the  fluted  condenser  tube  bundle  increases,  how¬ 
ever,  due  to  the  tube  enhancement.  This  may  be  off-set 
by  foundation  and  structural  benefits,  but  that  is  not 
addressed  here.  With  a  change  in  material.  Table  V-l 
shows  that  both  volume  and  weight  can  be  saved  at  the 
cost  of  increased  pumping  requirements. 

Heat  transfer  coefficients  for  condensation  are 
predicted  by  the  method  described  in  Chapter  III.  For 
flutes  which  are  of  the  size  12-24  mils  in  amplitude, 
increases  in  the  heat  transfer  coefficients  for  con¬ 
densation  are  greater  than  eight  times  that  predicted 


Table  V-2 


Comparison  Between  a  Horizontal  and  a  Fluted  Condenser* 


Horizontal 

Vertical 

Do  (in) 

.625 

.625 

t  (in) 

.065 

.065 

a  (mils ) 

- 

18 

a/p 

- 

.13869 

(°F) 

65.0 

65.0 

vc  .( ft/sec) 

8.0 

8.0 

w  (GPM) 
c 

7900 

7900 

No.  of  passes 

2 

2 

nt 

3294 

3294 

P_  (in-Hg-abs) 
c 

5.0 

5.0 

h__  (BTO/hr2oF) 
sc 

2000 

2000 

K  (BTU/hrft°F) 

17.0 

17.0 

L  (ft) 

14.85 

12.0 

H  (BTU/hrft2°F) 
o 

435.8 

538. 6(U  =509. 
n 

It  (°f) 

12.07 

1.639 

(q/A0)  (BTU/hr  ft2) 

19627 

24257 (q/An= 
2293C 

tc2  (°F) 

106.15 

106.15 

^total 

1.57xl08 

1.57xl08 

Bundle  Dia.  (in) 

53.4 

61.8 

Bundle  wt  (tons) 

9.63 

10.31 

•Material  for  both  is  Cu-Ni  70-30 
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by  the  Nusselt  analysis  for  horizontal  tubes.  Stress 
concentrations  nay  be  a  concern  for  submarine  use,  but 
by  limiting  the  amplitude-to-pitch  ratio,  stresses  can 
be  kept  below  an  acceptable  yield  criterion.  Great 
gains  in  weight  and  volume  savings  could  be  realized 
with  enhancement  on  the  cooling  water  side  of  the  tube 
in  addition  to  vertical  flutes. 

D.  Recommendations 

Because  thermal  resistance  is  treated  in  series, 
the  gains  made  on  the  condensation  side  of  a  tube  may 
not  be  fully  realized  because  of  the  insensitivity  of 
cooling  water  flow  and  resistance  due  to  scaling  or 
fouling . 

Further  investigation  is  required  to  consider  the 
following: 

1.  Internal  enhancement:  Tremendous  gains  may  be 
realized  by  increasing  the  heat  transfer  co¬ 
efficient  of  the  coolant.  An  internal  fin 
arrangement  could  accomplish  this.  Drawbacks 
to  such  arrangements  are  susceptibility  to 
fouling  and  tube  inspectability . 


2.  Bio fouling:  The  compactness  of  a  titanium 
marine  condenser,  or  possibly  an  aluminum 
stationary  condenser,  may  be  to  no  avail  if 
they  are  rapidly  blocked  or  fouled  by  marine 
organisms.  Active  and  passive  protection 
systems  must  be  investigated  which  retard  or 
prevent  biofouling,  and  if  these  are  net 
feasible,  an  in-service  cleaning  method  must 
be  considered. 

3.  Scaling:  This  allowance  is  a  bona  fide  guess. 
With  new  materials,  such  as  titanium,  investi¬ 
gation  into  scaling  allowances  and  mechanisms 
would  allow  a  less  conservative  approach  to 
condenser  design. 

4.  Experimental  Data:  This  is  an  absolute  neces¬ 
sity.  None  has  been  provided  with  this  thesis 
because  very  little  is  available  in  the  liter¬ 
ature  concerning  fluted  tubes.  A  criticism  of 
the  literature  cited  is  that  no  common  termi¬ 
nology  is  used,  nor  is  design  or  zero  point  data 
usually  provided.  Furthermore,  most  work  on 
fluted  tubes  has  been  performed  with  ammonia 
and  various  regrigerants . 


E .  Summary 

A  proposal  has  been  made  in  this  thesis  to  design 
a  lighter  more  compact  marine  condenser.  Calculations 
have  been  presented  for  selecting  different  condenser 
tubes  for  submarine  use.  Alternative  materials  have 
been  considered  and  a  comparison  condenser  design  has 
been  presented.  A  methodology  for  calculating  the  heat 
transer  coefficient  for  fluted  tubes  has  been  presented, 
and  this  is  incorporated  into  the  condenser  sizing  pro¬ 


cedure 


APPENDIX  A 


Condenser  Tube  Stress  Calculations 


Nomenclature 

a 

d 

h 

p 

p 

t 

Y 

YS 

a 

a 

r 

9  > 

Z 


amplitude  of  the  flute  (in) 
tube  diameter  (in) 
design  depth  (ft) 
pitch  of  the  flute  (in) 
pressure  (psi) 
tube  wall  thickness  (in) 
yield  criterion  (ksi) 
yield  stress  (ksi) 
stress  intensity  factor 
stress  (ksi) 

cylindrical  coordinates 


Subscripts 

i  inner  diameter 

o  outer  diameter 


x 


extreme  diameter 
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The  general  equations  for  stress  in  a  thick  walled 
aon-f luted  tube  are  £35] : 


-  *i[jva)2  -  i]  +  p0[(Vdi)2  -  (Vd)2. 

(Vdi)2  -  1 

E’i[(Vd)2  +  ’•j  -  p0[(Vdi)2  +  (do/d)2_ 

(<VdJ2  -  1 


a  a  -  P.  (this  is  because  the  tube  sheet  (A. 3) 

z  i 

distributes  external  pressure  on 
the  ends  of  the  tubes) . 

These  stresses  are  shown  in  Figure  A-l.  For  a  condenser 
PQ  <<  P^  and  equations  (A.l)  and  (A. 2)  simplify  to: 

-Pi[(Vd)2  -i] 

Or  *  -  (A. 4) 

(Vdi) 2  -  1 

piM2  +  d 

- - 

(d^d^2  - 1 


(A. 5) 


IX 


WB.TM5 


Trar’T*’: 


.  1  -  ~  '  *  v »'  ■  r- 
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Maximum  stresses  occur  at  the  inner  diameter  when  d  =  d^ 


max 


max 


Pi 


Pi 


1  + 


(<Vdi)  -  1 


(A. 6) 

(A.  7) 


a  *  -  P,  (A. 8) 

Z  l 

The  failure  criterion  used  is  the  Mises  yield 
criterion,  Y  [35,36]  where 

Y  *  {l  0r-«e)2  +  2  +  (az“ar)2]}  /  (A*9) 

For  a  smooth  surface  tube  this  reduces  to 

(A. 10) 


y  «  2  p. 


l  + 


(Vdi)  - 1 


or 


Pi  + 


3 


max 
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For  a  fluted  tube,  a  stress  concentration  occurs  at 
the  bottom  of  the  flute  valley.  The  mechanical  restriction 
to  selecting  a  fluted  tube  geometry  is  to  have  maximum 
stresses  in  the  fluted  tube  which  do  not  exceed  a  yield 
criterion  which  is  no  greater  than  the  yield  criterion  for 
a  smooth  tube.  Since  ur  =  0  at  the  surface  of  the  tube, 
the  Tresca  yield  criterion  is  appropriate  in  the  form; 


03 

° flute 
2 


(A. 11) 


or 


+  a, 


flute 


The  tangential  stress  for  the  flute  tube  at  the  bottom  of 
the  flute,  where  a  -  dQ,  is 


.Pi[(Vao)2  +1] 

(Vdi)2  -  1 


(A. 12) 


and  dx  =  dQ  +  4a 


(A. 13) 


a  is  a  stress  intensity  or  stress  concentration  factor.  3y 
requiring  the  yield  criterion  in  equation  (A. 11)  to  be  less 
than  or  equal  to  the  yield  criterion  in  equation  A-2, 
a  relationship  between  maximum  stress  for  the  two  tubes  is 


a 


9  flute 

0  9 

max 


<  1 


(A. 14) 


For  a  given  fluted  tube,  compared  to  a  smooth  tube  of  the 
same  inner  and  outer  diameters,  a  maximum  value  for  a  can 
be  determined.  Neuber  [37]  shows  a  method  to  relate  the 
amplitude-to-pitch  ratio,  a/p,  with  the  stress  intensity 
factor.  Thus  a  maximum  value  for  a/p  can  be  determined. 
This  concentration  factor  for  a  sine  flute  can  be  shown  to 
be 


(A. 15) 


and 


r 


o 


a 

( 2ira/p)  2 


(A. 16) 


Therefore , 


a  *  1  +  ^  2tt  /2 y 
P 


(A. 17) 
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Y  is  the  load  relieving  factor  which  is  described  in 
reference  [37]  and  accounts  for  the  presence  of  many  flutes 
on  the  surface  of  the  tube.  It  is  presented  in  Figure  (A-2) . 

Y  cam  be  related  to  some  mechanical  property  of  the 
tube  material  such  as  yield  stress  or  fatigue  limit.  For 
a  submarine,  a  very  conservative  safety  factor  of  5  can  be 
assumed  [38].  Therefore, 

Y  <  .2YS  (A. 10) 

From  reference  [2] 


*  h/2248  Ksi 

with  h,  the  design  depth  in  feet. 
Then  for  a  smooth  tube 


and  for  a  fluted  tube 


(A. 19) 


(A. 20) 
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Both  of  these  equations  contain  the  constraints  due  to 
material  strength,  design  depth,  and  tube  geometry.  The 
tube  wall  thickness  for  a  smooth  tube  can  be  readily  cal¬ 
culated  from  equation  (A. 20),  and  the  flute  limitations 
can  then  be  determined  from  equations  (A. 19)  and  (A. 17) 
which  will  ensure  that  the  yield  criterion  is  not  violated 
due  to  stress  concentrations. 

An  example  calculation  is  provided. 


EXAMPLE 

Two  different  materials  are  considered  for  a  2000  ft 
design  depth  condenser.  These  are  Cu-Ni  70-30  and  Ti. 
Other  specifications  are: 


do 

S 

.625  in 

a 

- 

.018  in 

d  ' 

X 

- 

.697  in 

YS (Cu-Ni) 

- 

25  Ksi 

YS (Ti) 

x 

40  Ksi 

From  equation  (A. 20) ,  the  inner  diameter,  d^,  can  be 
calculated  and  hence  the  wall  thickness,  t.  Wall  thick¬ 
ness,  however,  is  rounded  off  to  the  nearest  even  gauge. 


wks  *tsj*a  o  cr  ct 


V 
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.  l  A  A  l"l  *.' 


Cu-Ni 

Ti 

ai 

(in) 

.5016 

.5511 

t 

(in) 

.0617 

.0369 

di 

(rounded  off) 

(in) 

.495 

.555 

t 

(rounded  off) 

(in) 

.065 

.035 

Using  the  rounded  values  in  equation  (A. 21)  to  calculate, 
a,  and  then  using  a  to  determine  maximum  a/p  gives 


a 

2.0235 

2.0559 

y 

CM 

00 

• 

.81 

a/p 

.13 

.13 

In  light  of  the  conservative  allowance  for  yield  criterion 
select  a/p  in  both  cases  as  a/p  *  .15. 
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APPENDIX  3 

Calculations  for  Condensate  Film  Thickness 

Nomenclature 

a  amplitude  of  the  flute  (ft) 

A  area  (ft2) 

B  dimensional  group  defined  in  equation  (B.6) 

.  .  2 
gc  gravitational  cons tent  (lbm  ft/lbf  hr  ) 

h  heat  transfer  coefficient  (BTU/hrft2°F) 

hfg,  latent  heat  (BTU/lbm) 

K  thermal  conductivity  (BTU/hrft°F) 

Nu  Nusselt  number 

p  pitch  of  the  flute  (ft) 

P  pressure  (lbf/ft2) 

q  heat  flow  (BTU/hr) 

r  radius  of  curvature  of  the  fluted  sur¬ 

face  (ft) 

Sv  flute  half  perimeter  (ft) 

Sc  flute  arc  length  on  which  condensing 

occurs  (ft) 

T  temperature  (°F) 

AT  Ts  -  Tw  (°F) 

velocity  (ft/hr) 


V 


tube  surface  coordinates 


Greek  letters 


Subscripts 


height  of  the  condensate  in  the  center 
of  the  flute  (ft) 

mass  flow  in  S  direction  per  unit  length 
(lbm/hrft) 

condensate  film  thickness  (ft) 
angular  representation  of  a  point 
along  the  flute  surface 
dynamic  viscosity  (lbm/hrft) 
density  (lbm/ft3) 
surface  tension  (lbf/ft) 

non-dimensional  group  defined  in  ecn  (B.12) 
non-dimensional  group  defined  in  ecn  (B.6) 


coolant 

at  flooding  conditions 
inside 


nominal 


outside/initial  conditions 


luuuiu 


,,,  ■,,,  .in  ..  n  V'fi^PT 


R 


w 


Symbols 


T* 
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saturation 

wall 


average 


The  surface  considered  is  shown  in  Figure  III-2. 

The  following  assumptions  are  stated  explicitly  and  apply 
to  Figure  III-2, [7,9,12,13] . 

1.  The  flow  of  the  condensate  from  peak  to  valley 
on  the  flute  surface  is  thin  laminar  flow. 

2.  Because  the  film  is  thin,  intertia  forces  are 
negligible. 

3.  There  is  no  inter facial  shear  stress  between 
the  liquid  and  the  vapor. 

4.  There  is  no  flow  in  the  y  direction.  Flow  in 
the  z  direction  occurs  only  within  the  trough 
of  height  o. 

From  a  momentum  balance  on  the  surface  specified. 


p 


dP  _d_  fol 
dS  “  dS  (Fj 


Integrate  this  with  boundary  conditions 


(B.l) 


dVs 

-gj.  -  °  at  y  -  S 


and 


Vg  *  0  at  y  =  0 


The  equation  for  velocity  is 


V, 


1  dP 

u  3S 


4 


-  5y 


(B . 2a) 


and 


V  =  — 
S  5 


vSdy  " 


1_  dP  5__ 
u  dS  3 


(B.2b) 


Define  the  mass  flow  rate  per  unit  length. 


r  *  pvs5 


£_  dP  £2 

~  dS  3 


(B. 3a) 


and 


dr 

as 


d 

as 


£  dP  ^f] 
u  dS  3 


J 


(B. 3b) 


From  an  energy  balance 


q/A  -  m  -  hf  g 
s  o  fg  dS 


(3.4) 


dr 


Solve  eqn  (B.4)  for  ^5  and  set  equal  to  eqn  (B.3b) 


AD-A124  062  MARINE  SURFACE  CONDENSER  DESIGN  USING  VERTICAL  TUBES  2/2 

WHICH  ARE  ENHANOED(U)  NAVAL  POSTGRADUATE  SCHOOL 
MONTEREV  CA  C  G  BARNES  JUN  81 


UNCLASSIFIED 


F/G  12/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


■*.  r; 
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KAT 


fg 


d 

aS 


"  3  ,  n 

p  _d_  M  j 

u  3  dS  [r]  j 


which,  for  constant  physical  properties  becomes 


uKAT 

0<,h£g3c 


i  A 

3  dS 


63  <L 

dS 


Define  £2, 


a  2 


uKAT 


pchfggca 


BAT 

a 


So  eqn.  (B.5b)  may  be  written  as 


o  5  d 

Qa’33  3s 


,3  d  11 
5  dS  [r 


For  a  sine  function  flute  the  radius  of  curvature, 

2  ,  T3/2 


(*  *  (")•“’•)■ 
W 


cose 


The  arc  length,  s,  is 


» 

md 

EL  {i  +  I 

r2*al 

si 

4 

ST  \l  +  I 

i  p ) 

2  .  2  ^2  . 

sm  %  r  d<J> 


(3.5a) 

(B.5b) 

(B.  6) 

(B.5c) 

r ,  is 

(3.7) 

(B. 8a) 


as  shown  in  Figure  III-2. 


e 


c 


at  S. 


(B. 8b) 


For  given  values  of  a,  p,  and  0,  eqn.  (B.5c)  can  be  solved 
numerically  [7,12].  By  symmetry 


d  5 

as 


0  at 


3*0 


Eqn.  (B. 5c)  becomes 


S=0 

9*0 


3an 


(B .  9) 


which  gives  the  initial  value  for  the  boundary  layer 
thickness,  5Q. 


(B.10) 


Integrating  equation  (B.5C)  gives 


T'.-vf.-srT.vt^TV; v  «;■  V V* i".  n . '.T.  l»  y ? 


-yt>* 


,rs 


d 


[U 

IrJ 


3a.°. 


0 

f  4 


S  1 
dS  ( 
6 


dS 


Define  4  such  that 


and  by  finite  difference 


n+ 1 


+  4 


n 


where 


*Sn 


n+1 


S 


n 


Eqn.  (3.11)  may  be  written  in  the  form  of  a 
integration  to  give 


(3.11) 

(3.12a) 

(B. 12b) 

(B. 13) 

trapezoid  rule 

(B.  14a) 


which  becomes 


(3.14b) 


1  1 
*n  "  ?n+l 


3an 


2$  + 
n 


5n  cam  now  be  determined  for  A3  increments,  with  the  initial 
conditions 

&Q  from  eqn  (B.10) 

a  »  0 

o 

so  -  o 


Addendum  1  to  this  appendix  is  a  program  written 
for  a  programmable  calculator  to  solve  eqn.  (B.14b) 

The  average  film  thickness  of  the  condensate  can 
be  calculated 


S 


0 


Define  the  heat  tramsfer  coefficient,  h 


(B. 15) 


h 


(B. 16) 


The  average  heat  transfer  coefficient  is 


(B. 17a) 


AV.  A 


98- 


M 


,v*. 


where  Sc  is  shown  in  Figure  III-2  and  is  the  arc  length 
over  which  condensation  occurs.  By  definition 


K  *  . 1 31 , 

ns  isc<iT 


(B.  17b) 


where  q'  is  the  heat  transfer  per  unit  length  over  the 
entire  flute  perimeter.  The  heat  transfer  coefficient  on 
the  nominal  diameter,  hn,  may  be  obtained  from  the  equation 
for  heat  transfer 


q  *  h  A  A  T 


From  energy  conservation. 


(B. 18a) 


and  therefore 


Ha  «  Ha 
n  n  ss 


Then 


(B . 18b) 


(B. 18c) 


2S_ 

hn  *  “T  Hs  -  2K*<Sc> 


(B.  19) 


A.  Nusselt  number  may  be  defined  as 


Knp 

Nu  2  -g-  -  2«{SC) 


(B. 20a) 


and  for  S„  *  S„ 

C  V 


Nu  **  Nu  (S  )  *  2$(SV) 
o  v 


(B . 20b) 


From  equations  (B.6),  (B.8a),  (B.8b) ,  and  (III-10) 


Hn  -  f  (fl,  a/p,  w/wf) 


(B .  21) 


.*>  .  •  .*«  *.  .V 
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APPENDIX  B  -  ADDENDUM  1 
Calculator  Program  for  the 
Solution  of  Equation  (B.14b) 

This  program  was  written  for  a  Hewlett  Packard 
HP-41C  programmable  calculator.  It  requires  two  addi¬ 
tional  memory  modules,  and  a  printer  is  highly  desireable. 
Twenty  six  (26)  memories  must  be  allocated.  The  program 
is  run  by  simply  executing  the  program  label,  "DELTA". 

The  calculator  will  respond  by  inter-actively  asking  for 
inputs.  The  required  inputs  are: 

a,  p,  and  fl. 

The  outputs  are  for  n  »  0  to  n  *  60, 

r  ,  S  »  AS  .  5  ,  and  $  . 
n'  n'  sn'  n'  n 

The  method  of.  solution  is  exactly  as  explained  in 
Appendix  B.  Equation  (B.8a)  is  solved  by  a  10  point 
Gauss -Legendre  quadrature  method  for  a  sine  geometry  flute. 


i 


r  IN'.'. 


77  n  *;a  ;*  »» -  .  v 


*T7  ryrry  r*v~y 
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01 

LBL  *  DELTA1 * 

02 

"COMPUTER  BOUNDA1 

03 

"DRY  LAYER" 

04 

AVXEW 

05 

"AS  A  FCT.  OF  AR 

06 

"C  LENGTH . " 

07 

AVIEW 

08 

ADV 

09 

I************** 

10 

AVIEW 

11 

"S<N>/PHI<N> 

-DEL 

12 

"TA<tT>,  AVG1’ 

13 

AVIEW 

14 

••************<• 

15 

AVIEW 

16 

ADV 

17 

CLRG 

18 

SCI  4 

19 

1.488743390 

2-1 

20 

STO  10 

21 

2.955252247 

E-l 

22 

STO  11 

23 

4.333953941 

E-l 

24 

STO  12 

25 

2.692667193 

E-l 

26 

STO  13 

27 

6.794095683 

E-l 

28 

STO  14 

29 

2.190863625 

E-l 

30 

STO  15 

31 

8.650633667 

E-l 

32 

STO  16 

33 

1.494513492 

E-l 

34 

STO  17 

35 

9.739065285 

E-l 

36 

STO  18 

37 

6.667134431 

E-2 

38 

STO  19 

39 

GTO  10 

Clear  all  memories 


Loading  into  memory 
locations  weighting 
values  necessary  for 
integration  by 
quadrature 


*.v> 
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40 

LBL 

11 

41 

"INTEGRATIO 

42 

STO 

21 

43 

XOY 

44 

STO 

22 

45 

- 

46 

2 

47 

/ 

48 

STO 

22 

49 

RCL 

22 

50 

RCL 

21 

51 

+ 

52 

2 

53 

/ 

54 

STO 

21 

55 

0 

56 

STO 

00 

57 

10.021 

58 

STO 

25 

59 

XEQ 

12 

60 

XEQ 

12 

61 

XEQ 

12 

62 

XEQ 

12 

63 

XEQ 

12 

64 

RCL 

00 

65 

RCL 

20 

66 

* 

67 

RTN 

68 

LBL 

12 

69 

"QUADRATURE 

70 

"ROUTINE” 

71 

RCL 

IND 

25 

72 

ISG 

25 

73 

STO 

22 

74 

CHS 

75 

RCL 

20 

76 

* 

77 

RCL 

21 

78 

+ 

79 

XEQ 

14 

80 

RCL 

IND 

25 

81 

* 

82 

ST+ 

00 

83 

RCL 

22 

84 

RCL 

20 

85 

* 

86 

RCL 

21 

87 

+ 

88 

XEQ 

14 

89 

RCL 

IND 

25 

90 

* 

91 

ST+ 

00 

92 

ISG 

25 

93 

RTN 

Integration  by  Gauss- 
Legendre  Quadrature 


b 

* 

f(x)d  * 
J  x 

a 


b-a 

~T~ 


10 

l  wif 

i-1  x 


(b-a) +b+a) ) 
- 5 - 

*  j 


I 

31 


94  LBL  14 

95  "INTEGRAND" 

96  SIN 

97  RCL  23 

98  * 

99  2 
100  * 

101  PI 

102  * 

103  RCL  24 

104  / 

105  X+2 

106  1 

107  + 

108  SORT 

109  RCL  24 

110  * 

111  2 
112  / 

113  PI 

114  / 

115  RTN 


Integrand : 


2-  (±  +  flial ' 
l*  \  l  P 


116  LBL  10 

117  "INITIALIZATION" 

118  "ROUTINE" 

119  RAD 

120  "A<FT>  ?" 

121  PROMPT 

122  STO  23 

123  "P<PT>?" 

124  PROMPT 

125  STO  24 

126  "OMEGA?" 

127  PROMPT 


Data  input:  a 


3 

* 


Calculation  of  initial 
values  (continued) 


RCL 

23 

* 

STO 

09 

PI 

60 

/ 

STO 

08 

0 

STO 

02 

STO 

05 

STO 

06 

STO 

26 

RCL 

23 

RCL 

24 

/ 

2 

* 

PI 

* 

X+2 

1/X 

RCL  23 
* 

STO  04 
"RO<FT>* " 
ARCL  X 
AVIEW 
0 

"SO<FT>*" 
ARCL  X 
AVIEW 
FIX  3 
"PHIO-" 
ARCL  X 
AVIEW 


LB  =  tt/60 


rQ  from  equation  (B.7) 

9*0 

so  -  0 


Calculation  of  initial 
values  (continued) 


165  RCL  23 

166  RCL  24 

167  / 

168  2 

169  * 

170  PI 

171  * 

172  X+2 

173  3 

174  * 

175  1 

176  + 

177  RCL  23 

178  * 

179  1/X 

180  RCL  09 

181  * 

182  SQRT 

183  SORT 

184  RCL  24 

185  * 

186  2 

187  / 

188  PI 

189  / 

190  STO  03 

191  SCI  4 

192  "DELTAO<FT>='"  5  from  equation  (B.1C) 

193  ARCL  X  ° 

194  AVIEW 

195  2 


198  LBL  01 

199  A DV 

200  DSE  25 

201  RCLV  25 

202  STO  07 

203  1 

204 

205  SF  12 

206  FIX  0 

207  "N»" 

208  ARCL  X 

209  AVIEW 

210  CF  12 

211  SCI  4 

212  XEQ  03 

213  STO  26 

214  STO  04 

215  "R<N>  <FT>=" 

216  ARCL  X 

217  AVIEW 

218  RCL  26 

219  RCL  04 

220  X<>Y 

221  STO  04 

222  1/X 

223  CHS. 

224  X<>Y 

225  1/X 

226  + 

227  STO  01 

228  0 

229  ENTER+ 

230  RCL  07 

231  ENTER* 

232  RCL  25 

233  STO  07 

234  RON 

235  XEQ  11 

236  STO  27 

237  RCL  05  s 

238  "S< N><FT>="  n 

239  ARCL  X 

240  AVIEW 

241  RCL  27 

242  RCL  05 

243  X<>  Y 

244  STO  05 

245 

246  CHS 

247  "DEL/ S<N>  <FT>=” 

248  ARCL  X 

249  AVIEW 


Solution  to  equation  (B.14b) 


n 


Subroutine  for  equation  (B.7) 


rn  from  equation  (B.7) 


Integration  subroutine 

nAS 

*  [  f (9)de  from  equa- 

/.  tion  (3.8a) 


ASn  from  equation  (3.13) 
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250  STO  06 

251  RCL  07 

252  STO  25 

253  RCL  03 

254  STO  07 

255  SCI  4 

256  FC?C  01 

257  XEQ  02  Subroutine  to  iterate  for 

258  "DELTA <N>  <FT>*"  5  from  ecuation  (B.14b) 

259  ARCL  x  n 

260  AVIEW 

261  RCL  06 

262  RCL  07 

263  / 

264  RCL  02  $  from  ecuation  (B.12b) 

265  FIX  3  n 

266  "PHI<N>=" 

267  ARCL  X 

268  AVIEW 

269  + 

270  STO  02 

271  2 

272  ST+  25 

273  RCL  25 

274  63 

275  X-Y? 

276  GTO  05 

277  GTO  81 
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312  LBL  03 

313  RCL  08 

314  ST*  07 

315  RCL  07 

216  SIN 

317  RCL  23 

318  2 

319  * 

320  PI 

321  * 

322  RCL  24 

323  / 

324  STO  01 

325  * 

326  Xt2 

327  1 

328  + 

329  1.5 

330  YtX 

331  RCL  07 

332  COS 

333  / 

334  RCL  01 

335  X+2 

336  / 

337  RCL  23 

338  * 

339  RTN 

340  LBL  05 

341  BEEP 

324  ADV 

343  ************** 

344  AVIEW 

345  ADV 

346  RTN 

347  STOP 

348  .END. 


Subroutine  for  eauation 
(B .  7) 

rn  from  equation  (B.7) 
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AP  PEND IX  c 

Fluted  Condenser  Sizinc  Procedure 
■  — 1  . 


Nomenclature 

a  amplitude  of  the  flute  (ft) 

A  area  (ft2) 

B  dimensional  group  defined  in  Appendix  B 

D  diameter  (ft) 

F  number  of  flutes  on  the  tube  surface 

2 

g  gravitational  constant  (lbm  ft/lbf  hr  ) 

C 

h  heat  transfer  coefficient  (BTU/hrft2°F) 

latent  heat  (BTU/lbm) 
k  thermal  conductivity  (BTO’/hrft°F) 

1  ligament  (ft) 

L  length  (ft) 

M  non-dimensional  group  defined  in  Chapter  III 

N  number  of  tubes 

Nu  Nusselt  number 


p  pitch  of  the  flute  (ft) 

P  pressure  (in-Kg-abs) 

q  heat  flow  (BTU/hr) 

t  tube  wall  thickness  (ft) 

T  temperature  (°?) 

AT  T-T  ( °F) 

3  W 

U  overall  heat  transfer  coefficient 

(BTU/hrft2°F) 


V  velocity  (ft/hr) 

w  axial  mass  flow  of  liquid  (lbm/hr) 

Wt  weight  (lbm  or  tons) 

Greek  letters 

u  dynamic  viscosity  (lbm/hrft) 

p  density  (lbm/ft3) 

a  surface  tension  (lbf/ft) 

8  non-dimensional  group  defined  in  Chapter  III 

Subscripts 

B  tube  bundle 

c  coolant 

D  tube  bundle  diameter 

f  at  flooding  conditions 

i  inside 

n  nominal 

o  outside 

s  saturation 

T  total 

w  wall 

sc  scale 

x  cross  section 

Symbols 


average 


This  procedure  will  consist  of  an  example  which 
sizes  a  submarine  condenser.  As  required  in  reference 
[11] ,  a  submarine  condenser  has  an  even  number  of  passes 
for  the  coolant  flow.  A  maximum  length  will  be  specified 


Design  Specificatons 
Design  depth 
Material 


Tube: 


Coolant: 


o 

t 

a 

a/p 

F 

Dn 

Sea  water 
Tc, 


Wc 

Mo.  of  passes 
No.  of  passes 
Condensate : 

Pressure 

T 

1sat 

p  . 


2000  ft 

Cu-Ni  70-30 

5/8  in  (.625  in) 

.065  in  (16  Gauge) 

18  mils  (1.50  x  10-3  ft) 

.31205 

36 

.643  in 

65 . 0°F 
8.0  ft/sec 

7900  GPM  (4.0559  X  106  lbm/hr) 
2 

3294 


5.0  in-Hg-abs 
133. 7°F 

61.554  lbm/ft3 


.376  BTU/hrft°F 
.00459  lbf/ft 
1017.4  BTU/lbm 


u4  *  1.20  lbm/hrft 

Thermal  conductivity  of  tube  wall  *  17.0  BTU/hrft°F 
hscale  “  2000  BTU/hrft2°F 

Total  tube  length  per  pass  *  12.0  ft 
Preliminary  calculations: 

Nuon1(/4  ■  4.612  from  Figure  III-4 

B  -  3.7649  x  10“12  ft/#F  eqn.  (III-7) 

Wf  -  6.7157  lbm/hr  eqn.  (III-10) 

The  condenser  will  be  segmented  into  6  sections  as 
shown  in  Figure  (IV-1) .  The  length,  L,  of  section  1 
will  be  that  length  required  for  condenstate  run-off  of 
w/wf  ■  0.060.  This  length  will  be  used  for  sections  2, 
5,  and  6.  The  lengths  for  sections  3  and  4  will  be 
the  length  necessary  to  make  each  pass  total  length 
equal  to  12.0  ft. 


A; _ Calculation  of 

1.  Specify  w/wf 

then  M  *  .0820  from  equation  (III-19) 

2.  Assume  TC2  *  66°F;  i  *  1 

3.  Assume  L  »  10.0  ft;  i  1 

a.  Solve  for  from  equation  (III-29) 

b.  Solve  for  (q/An)^  from  equation  (HI-25) 

c.  Solve  for  from  equation  (HI-23) 

d.  Repeat  a,  b,  and  c  until  +  ^ 

4.  Solve  for  Tc2j+i  from  equation  (III-24) 

5.  Repeat  3  and  4  until  Tc2j  *  Tc2j+i 
Results  are  shown  in  Table  C-l. 


B.  Calculation  of  (w/wf]_2 

1.  Specify  L2  »  Lx  -  4.957  ft. 

2.  Assume  Tc2  ■  77.36°F;  j  ■  1 

3.  Assume  AT  *  0.50°F;  i  *  1 

a.  Solve  for  C<?/An) ^  from  equation  (III-25) 

b.  Solve  for  (w/w^)^  from  equation  (III-23) 

c.  Solve  for  M  from  equation  (III-19) 

d.  Solve  for  ATi+^  from  equation  (III-29) 


e.  Repeat  a-d  until  at  .  ■  AT. .. 


TABLE  C- 


MM  WMT* 


B u1.  »*.  » i  •  "*N  S’*  '.*7 
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4.  Solve  for  Tc2j+i  from  equation  (Hi-22) 

5.  Repeat  3  and  4  until  Tc2j  =  Tc2j+i 
Results  sure  shown  in  Table  C-  2. 


Calculation  of  all  remaining  length  sections  pro¬ 
ceeds  as  for  section  2,  with  the  desired  length  speci¬ 
fied.  The  coolant  exit  temperature  for  each  section  is 
the  inlet  temperature  for  the  succeeding  section. 

Results  for  the  entire  condenser  are  shown  in  Table  C-3 . 

A  check  on  the  solution  to  equation  (III-22) 
verify  the  results  in  Table  C-3. 


E 


.6027 


hfq  wf 

a,2307[f (a/p)] *9226 

L  AT 

>  , 

P 

UAT 

2307 


(III-22) 

where  f[a/p]  «  Nuq01//4 
This  should  be  compared  to 
(q/A  ) 

— — —  *  E  shown  in  Table  C-3. 

AT 


Table  C-4  shows  this  comparison. 
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TABLE  C-4 

Comparison  Between  Equation  (III- 22] 
and  Calculated  Results 


Section 


hn  by  eqn  (22) 


h  by  iterative 
solution 


12322 


12312 


13206 


13197 


15313 


15302 


15812 


15800 


15190 


15178 


16452 


16441 


C.  Sizing  the  Tube  Bundle  Diameter 


The  number  of  tubes  across  the  diameter,  N^,  of 
a  circular  bundle  containing  total  number  of  tubes  is 


(C.l) 


The  geometric  diameter  for  the  individual  tube  used  in 
sizing  the  tube  bundle  is  the  extreme  diameter,  D.. 


Dx  -  0o  +  4a 


(C.  2) 


The  physical  separation  in  the  tube  sheet  between  outer 
edges  of  adjacent  tubes  is  termed  the  ligament,  1. 
Minimum  ligament  distance  is  0.212  in.  [42].  This  limit: 
is  set  by  manufacturing  constraints.  Thus  the  tube 
bundle  diameter,  D0,  is 

°B  “  ND  •  °x  +  (ND  -  1)  •  1  (C-3) 

For  the  condenser  considered, 

Nd  »  65 

Dx  *  .6970  in. 

Dg  *  58.87  in. 

Allowing  5  percent  for  steam  flow  passages, 
structural  members,  air  ejectors,  and  downcomers  gives 
an  external  diameter  for  the  condenser  bundle 


L2 


D.  Weight  Estimation  of  the  Tube  Bundle 

The  weight  of  an  individual  tube  can  be  deter- 


mined  by 

section. 

> 

the  tube 

WtB  -  Ax  L  p  Nt 


(C .  4 ) 


The  cross  sectional  area  for  a  sine  geometry  flute  is 
approximated  by 


(D2  -  D?) 
'  n  x' 


(C.  5) 


For  the  condenser  considered, 


Ax  »  .1507  in2 


WtB  -  23091  »  10.31  tons 


REFERENCES 


1.  Taggart,  R. ,  Ed.,  Ship  Design  and  Construction, 
Society  of  Naval  Architects  and  24arine  Engineers, 

New  York,  1980. 

2.  Jackson,  H.,  "Submarine  Design  Notes",  Dept.  XIII-A 
Professional  Summer  at  MIT,  1980. 

3.  Lecture  by  Mr.  Dan  Weiler  of  Naval  Sea  Systems  Com¬ 
mand,  3211,  to  MIT  XIII-A  students  on  8  January  1980. 

4.  "Shaping  the  General  Purpose  Navy  of  the  Eighties: 
Issues  for  Fiscal  Years  1980.-1985",  Congressional 
Budget  Office,  Congress  of  the  United  States, 

January  1980. 

5.  Couhat,  J. ,  Combat  Fleets  of  the  World,  1981/81 , 
United  States  Naval  Institute,  Annapolis,  MD,  1980. 

6.  Rohsenow,  W.  and  Choi,  H. ,  Heat ,  Mass  and  Momentum 
Transfer,  Prentice-Hall,  Englewood  Cliffs,  NJ,  1961. 

7.  Gregorig,  R. ,  " Hautkondensation  an  feingewellten 
Oberfl&chen  bei  Berucksichtigung  der  Oberflachen- 
spannung,"  Zeitschrift  fur  angewande  Mathematik 
und  Physik,  Vol.  5.,  1954,  pp.  36-49. 

8 .  Zener,  C.  and  Levi,  A.,  "Drainage  Systems  for  Con¬ 
densation,"  Journal  of  Engineering  for  Power, 

Vol.  96,  July,  19^4 ,  pp.  209-215.' 

9.  Webb,  R. ,  "A  Generalized  Procedure  for  the  Design 
and  Optimization  of  Fluted  Gregorig  Condensing 
Surfaces,"  Journal  of  Heat  Transfer,  Vol.  101, 

May,  1979,  pp.  335-3^9. 

10.  Conversation  with  Professor  B.  Mikic,  MIT  on 
26  February  1981. 

11.  Military  Specification,  "Condensers,  Steam,  Naval 
Shipboard,"  MIL-C-15430J  (Ships),  19  June  1974. 

Yamamoto,  H.  and  Ishibachi,  T. ,  "Calculation  of 
Condensation  Heat  Transfer  Coefficients  of  Fluted 
Tubes,"  Heat  Transfer  Japanese  Research,  Vol .  6 . , 
October  -  December,  1977*,  pp.  61-68. 


12. 


13 


Panchal,  C.  and  Bell,  K. ,  "Analysis  of  Nusselt- 
Type  Condensation  on  a  Vertical  Fluted  Surface," 

.  Condensation  Heat  Transfer,  ASME,  New  York,  1979, 
pp.  45-54 . 

14.  Mori,  Y.,  Hyikata,  K. ,  Hirasawa,  S.,  and  Nakayama, 

W. ,  "Optimized  Performance  of  Condensers  with  Out¬ 
side  Condensing  Surface,"  Condensation  Heat  Transfer 
ASME,  New  York,  1979,  pp.  55-62. 

15.  Krassik,  I.,  Krutzsch,  W. ,  Fraser,  W. ,  and  Messina, 
J.,  eds..  Pump  Handbook,  McGraw-Hill,  New  York,  1976 

16.  Harrington,  R. ,  ed..  Marine  Engineering,  SNAME, 

New  York,  1971. 

17.  Letter  from  Commander,  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  to  Commander,  Naval 
Sea  Systems  Command  (SEA  05DC3) ,  Ser  10310,  dated 

7  January  1980. 

18.  Metals  Handbook,  Ninth  Edition,  Vols.  2  and  3, 
American  Society  Ifor  Metals,  Metals  Park,  OH,  1979. 

19.  Efird,  K.,  and  Anderson,  D.,  "Sea  Water  Corrosion 
of  90-10  and  70-30  Cu-Ni:  14  Year  Exposures," 
Materials  Performance,  Vol.  14,  November,  1975, 
pp.  37-40. 

20.  Bulow,  C. ,  "Use  of  Copper  Base  Alloys  in  Marine 
Services,"  Naval  Engineers  Journal,  Vol.  77,  June 
1965,  pp.  476-432. 

21.  Danek,  G. ,  "The  Effect  of  Sea  Water  Velocity  on 
the  Corrosion  Behavior  of  Metals,"  Naval  Engineers 
Journal,  Vol.  78,  October,  1966,  pp.  763-769. 

22.  Fontana,  M.  and  Greene,  N. ,  Corrosion  Engineering, 

2d  ed. ,  McGraw-Hill,  New  York,  1978. 

23.  Litvin,  D.,  and  Smith,  D.,  "Titanium  for  Marine 
Applications,"  Naval  Engineers  Journal,  Vol.  83, 
October,  1971,  pp.  37-43. 

24.  Knudsen,  J. ,  "Apparatus  and  Techniques  for  Measure¬ 
ment  of  Fouling  of  Heat  Transfer  Surfaces,"  Con¬ 
denser  Biofouling  Control,  Symposium  Proceedings, 
Electric  Power  Research  Institute,  Ann  Arbor 
Science  Publishers,  Ann  Arbor,  MI  1980. 


-124- 


Characklis,  W. ,  Bryers,  J.,  Trulear,  M.  and  Zelver, 
N.,  "Biofouling  Film  Development  and  Its  Effects  on 
Energy  Losses:  A  Laboratory  Study,"  Condensor  Bio¬ 
fouling  Control,  Symposium  Proceedings,  Electric 
Povfsir  Research  Institute,  Ann  Arbor  Science  Pub¬ 
lishers,  Ann  Arbor,  MI,  1980. 

Letter  from  Commander,  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  to  Commander,  Naval 
Sea  Systems  Command  (SEA  532),  Ser  2722,  dated 
23  July  1980. 

Bell,  K. ,  "The  Effect  of  Fouling  on  OTEC  Heat  Ex¬ 
changer  Design,  Construction,  and  Operation,"  Un¬ 
published,  received  by  private  communication, 
September  1980. 

Blume,  W.  and  Kirk,  B. ,  "Application  of  the  Cathelco 
Antifouling  System  for  the  Control  of  Marine  Growth 
in  Seawater  Inlets  and  Condensers , "  Condenser  Bio¬ 
fouling  Control,  Symposium  Proceedings,  Electric 
i>ower  Research  Institute,  Ann  Arbor  Science  Pub¬ 
lishers  ,  Ann  Arbor,  MI,  1980. 

Conversation  with  Mr.  Alan  Fritsch,  David  W.  Taylor 
Naval  Ship  Research  Development  Center,  Annapolis, 
Maryland,  on  14  June  1980. 

Brown,  B.,  et  al.  Marine  Corrosion  Studies,  NRL 
Memorandum  1574,  Second  Interim  Report  of  Progress, 
U.S.  Naval  Research  Lab,  Washington,  DC,  November, 
1964. 

Uhlig,  H.,  Corrosion  and  Corrosion  Control,  2d  ed. 
John  Wiley  and  Sons,  New  York,  1971. 

Cotton,  J.  and  Downing,  B.,  "Resistance  of  Titanium 
to  Sea  Water,"  Transactions  of  the  Institute  of 
Marine  Engineers,  Vol.  69,  August,  1957,  pp.  311-319 . 

LaQue,  F.,  "Deterioration  of  Metals  in  an  Ocean 
Environment,"  Ocean  Engineering,  Vol.  1,  March, 

1969,  pp.  299-im 

Conversation  with  Mr.  Kurt  Bredehorst  and 
Mr.  Walter  Aerni  of  NAVSEA  on  15  June  1980. 


-125- 


35.  Crandall,  S.,  Dahl,  N.  and  Lardner,  T. ,  eds.. 

An  Introduction  to  the  Mechanics  of  Solids,  2d  ed. , 
NcGraw-Hill ,  New  York-,  1372. 

36.  McClintock,  F.  and  Argon,  A.,  eds..  Mechanical 
Behavior  of  Materials,  Addison-Wesley ,  Reading , 
Massachusetts,  1$66. 

37.  Neuber,  H. ,  Theory  of  Notch  Stresses,  Principles 
for  Exact  Stress  Calculation,  j.  w. Edwards , 

Ann  Arbor,  Michigan,  1946 . 

38.  Conversation  with  CAPT  Harry  A.  Jackson,  USN  (Ret.) 
on  18  November  1980. 

39.  Welty,  J. ,  Engineering  Heat  Transfer,  John  Wiley 
and  Sons,  New  York,  1§7S. 

40.  Fugii,  T.  and  Honda,  H. ,  "Laminar  Filmwise  Con¬ 
densation  on  a  Vertical  Single  Fluted  Plate," 

Sixth  International  Heat  Transfer  Conference,  Vol .  2 , 
National  Research Council  of Canada,  Toronto, 

Canada,  August  7-11,  1978,  pp.  419-424. 

41.  Harrington,  R. ,  ed..  Marine  Engineering,  Society  of 
Naval  Architects  and  Marine  Engineers ,  Mew  York, 

1971. 

42.  Conversation  with  Mr.  Henry  Brawn  of  DeLeval  Corp. , 
Florence,  NJ,  on  12  August  1980. 

43.  Titus,  M. ,  "Investigation  of  Rolled  Titanium  Tube- 
To-Tube  Sheet  Joints  for  Condenser  Application," 

David  W.  Taylor  Naval  Ship  Research  and  Development 
Center,  Bethesda,  MD,  October,  1978. 


END 

v 

a 

FILMED 

11-83 


